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predicting  properties,  as  is  possible  for  secondary  explosives. 

Essential  primary  explosive  properties  are  the  ability  to  remain 
ignited  and  quickly  undergo  deflagration  to  detonation  transition  (DDT). 
Several  practical  primaries  such  as  lead  styphnate  and  lead  azide  have 
detectable  burn  times  before  detonation.  Thus  the  study  of  DDT  in  cast 
explosives  was  initiated  to  establish  whether  DDT  ineasurements  were  use- 
ful as  a method  to  screen  explosives  for  primary  properties.  DDT  studies 
of  two  potential  castable  primaries,  the  isomers  l-methyl-5-nltrotetrazole 
(1-MNT)  and  2-methyl-5-nitrotetrazole  (2-MNT)  were  undertaken  as  a first 
step  to  achieving  this  objective. 

-^OT  experiments  showed  that  DDT  occurs  more  readily  in  cast  2-MNT  than 
in  cast  1-MNT. ^ The  threshold  voltages  for  ignition  by  hot  bridgewire  were 
found  to  be  110  volts  for  cast  2-MNT  and  375  volts  for  cast  1-MNT.  Critical 
times  and  distances  for  DDT  in  cast  2-MNT  of  8 /isec  and  2-3  mm,  respectively, 
; were  determined  by  performing  experiments  on  charges  of  different  lengths. 

' Critical  parameters  for  cast  1-MNT  were  not  determined  because  the  thresh- 
old voltage  for  Ignition  was  too  high. 

Macek's  model  for  DDT  in  cast  explosives  was  modified  in  a theoretical 
study  to  provide  a more  realistic  treatment  of  the  accelerating  flame. 
Calculations  were  made  with  the  deflagration  treated  as  a reactivity 
discontinuity  whose  burnt  state  satisfied  either  Chapman-Jouguet  (CJ) 
condition  or  the  condition  of  zero  particle  velocity  to  determine  upper 
and  lower  limits  for  the  deflagration  velocity  as  functions  of  pressure 
in  the  unburned  explosive.  >Times  and  distances  for  DDT  in  cast  pentollte 
and  crystalline  lead  azide  were  predicted  and  agreed  well  with  the  meager 
experimental  data  available. 

hm 

/ 


DDi 

EDITION  OF  1 NOV  6S  IS  OBSOLETE 


[(BACK) 


February  4,  1977 


Final  Report,  76-2 

December  5,  1975,  to  December  4,  1976 


structure/property  correlations  in  primary  explosives 

i 

By:  Craig  M.  Tarver,  Thomas  C.  Goodale, 

Michael  Cowperthwaite,  and  Marion  E.  Hill 


Prepared  for: 

U.S.  NAVAL  SEA  SYSTEMS  COMMAND 

Explosives  and  Pyrotechnics  Branch,  SEA  0332 

Washington,  O.C.  20300 

Attn:  Dr.  A.  B.  Amster 


Project  Serial  Number  SF3354-317 
NAVGSASYSCOM  Contract  Number  N00024-76-C5329 

SRI  Project  No.  PYU-4770 


Approved  by: 


'/ 1 


Vice  President 
ife  Sciences  Group 


CONTENTS 


FOREWORD tl 

LIST  OP  ILLUSTRATIONS  lit 

LIST  OF  TABLES  v 

I  INTRODUCTION  AND  SUMMARY  1 

II  STRUCTURAL  CONSIDERATIONS  6 

A.  Structure-Property  Characteristics  6 

B.  Structure-Property  Correlations 10 

1.  Metal  Ionization  Potentials  10 

2.  Bond  Rearrangement  11 

3.  Substltuency  Effects  In  Tetrazoles  14 

4.  Physical  Variables  16 

III  DEFLAGRATION  TO  DETONATION  IN  CAST  EXPLOSIVES  18 

A.  Background  18 

B.  DDT  Measurements  on  Cast  1-MNT  and  2-MNT 20 

C.  Modification  of  Macek's  Model  for  DDT  In 

Cast  Explosives  36 

D.  DDT  and  STD  Calculations  for  Lead  Azide 59 

IV  CONCLUSIONS  70 

REFERENCES  73 

APPENDICES 

A SYNTHESES  OF  1-  AND  2-METHYL-5-NITROTETRAZOLES  77 

B ESTIMATION  OF  HEATS  OF  FORMATION  AND  DENSITIES 

OF  1-  AND  2-METHYL-5-NITROTETRAZOLE  81 


1 


FOREWORD 


This  final  report  is  submitted  In  fulfillment  of  Contract  N00024- 
76-C-5329,  Navy  Project  Serial  Number  SF3354-316,  and  covers  the  period 
December  5,  1975  through  December  4,  1976.  For  completeness,  summaries 
and  excerpts  are  included  from  the  previous  two  contracts  on  the  same 
subject,  Contract  N00017-73-C-4346  and  Contract  N00024-75-C-5008. 

The  research  program  was  performed  by  staff  of  the  Physical  Sciences 
Division  under  the  supervision  of  Marlon  E.  Hill.  Craig  M.  Tar.ar  and 
Thomas  C.  Goodale  were  the  principal  experimentalists,  and  Michael 
Cowperthwalte  with  Craig  Tarver  developed  the  modified  theory  of  de- 
flagration to  detonation  transition. 

We  wish  to  acknowledge  the  extensive  and  helpful  discussions  and 
suggestions  by  our  monitors.  Dr.  A.  B.  Amster,  Head,  Explosives  and 
Pyrotechnics  Branch,  SEA  0332,  and  Mr.  Irving  Kabik  of  the  White  Oak 
Laboratory,  Naval  Surface  Weapons  Center. 
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I INTRODUCTION  AND  SUMMARY 


Primary  explosives  such  as  lead  azide  have  bemi  used  for  nearly  a 
hundred  years  to  Initiate  detonation  of  less  sensitive  secondary  explo- 
sives such  as  TNT,  Today  lead  azide  Is  the  principal  primary  explosive 
used  by  the  Department  of  Defense,  Because  of  several  deficiencies  with 
lead  azide,  notably  Its  unpredictable  sensitivity  often  arising  from 
copper  azide  formation,  many  other  energetic  materials  have  been  prepared 
for  the  same  purposes,  but  most  have  failed  to  find  general  acceptance  In 
the  explosives  community  for  one  reason  or  another. 

Consequently,  there  are  recurrent  requests  for  new  primaries  that 
are  not  In  the  azide  class.  Often  an  approach  to  fulfilling  the  need 
Is  to  search  the  shelves  for  any  explosive  that  one  suspects  Is  so  highly 
sensitive  that  It  may  be  a primary.  Such  an  approach  leads  inevitably 
to  an  Edlsonlan  try-and-see  series  of  syntheses  and  tests.  The  continuing 
use  of  lead  azide  Indicates  that  no  other  single  compound  has  been  found 
to  fulfill  all  requirements. 

Too  little  is  known  of  the  Influence  of  chemical  structure  in  a 
fundamental  way  on  the  action  of  initiators.  In  the  secondary  explosive 
field,  much  research  has  been  conducted  with  the  result  that  some  cor- 
relation between  chemical  structure  and  detonation  behavior  is  available 
to  assist  chemists  and  physicists  In  designing  and  applying  new  secondary 
systems.  Such  is  not  the  case  In  the  prlraary  field,  which  has  experimental 
results  and  some  empirical  correlations,  but  as  yet  no  general  theoretical 
basis. 

The  objective  of  this  research  was,  therefore,  to  place  the  choice 
of  primary  explosives  for  practical  applications  on  a firmer  theoretical 
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base  related  to  the  chemical  and  physical  properties  of  those  materials 
whose  action  Is  generally  classed  as  "primariness."  This  base  then  should 
permit  reasonable  correlation  between  structure  and  primary  properties, 
and  eventually  lead  to  the  design  of  explosive  molecules  suitable  for 
synthesis  and  testing  In  practical  detonation  compositions.  Replacement 
of  lead  azide  without  an  Edlsonlan  hlt-or-mlss  approach  would  be  feasible. 

Our  emphasis  was  placed  on  obtaining  as  much  information  as  pos- 
sible on  compounds  other  than  azides  since  this  class  has  had  rxtensive 
study  and  development.  (A  review  on  inorganic  azides  in  two  volumes 
Is  being  published  by  Plcatlnny  Arsenal.^)  Attempts  to  synthesize  new 
primary  explosives  have  concentrated  on  variations  of  metal  derivatives 
of  energetic  compounds,  for  example,  heavy  metal  compounds  such  as  lead 
styphnate  and  silver  azide  or  azo  compounds  such  as  dlazodlnltrophenol , 
tetracene  and  cyanurlc  triazlde.  Although  several  of  these  materials 
also  find  use  as  primaries,  their  applications  are  limited  and  they  are 
often  used  In  conjunction  with  lead  azide.  No  fully  organic  primaries 
are  used  at  this  time. 

The  very  special  requirements  for  efficient  primaries — sensitive 
but  not  too  sensitive,  fast  transition  from  deflagration  to  detonation, 
high  energy  yield,  and  reproducible  properties  when  prepared  in  production 
quantities— place  considerable  constraints  on  design  and  testing  of 
Initiators.  Guidelines  for  preparing  new  materials  to  meet  stringent 
requirements  for  useful  primary  explosives  are  based  largely  on  experience 
with  materials  that  work  rather  than  on  any  fundamental  principles  of 
kinetics,  thermochemistry,  and  physics.  We  believe  that  this  was  an 
opportune  time  to  attempt  vhe  characterization  of  primary  explosives  in 
a more  fundamental  way  than  had  been  possible  In  the  past  and,  on  the 
basis  of  these  relationships,  to  develop  new  superior  materials  and  ways 
of  testing  and  evaluating  them. 
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Kumsrous  research  workers  in  ^:he  field  of  explosives  were  Interviewed 
to  determine  the  state  of  the  art  and  to  Identify  the  major  unsolved 
problems  In  applications  of  primaries.  As  a complement  to  these  visits, 
we  examined  the  literature  for  gaps  In  knowledge.  Another  task  was  an 
attempt  to  correlate  the  physical  and  chemical  properties  of  explosives 
with  their  action  as  primaries.  The  results  were  not  quantitative  and 
provided  only  several  general  observations  and  conclusions  from  which 
an  approach  to  a possible  screening  method  was  derived. 

It  seems  much  easier  to  compare  a primary  explosive  with  a secondary 
explosive  than  to  define  a primary  explosive  In  absolute  terms.  Primary 
and  secondary  explosives  have  many  similarities.  Both  will  burn  to  deto- 
nation, although  secondaries  sometimes  go  out;  primaries  never  do.  Con- 
finement and  a "long"  burning  time  are  required  for  secondaries.  The 
existence  of  "near  primaries"  such  as  some  of  the  hydroxynltroaromatic 
compounds  and  "sensitive"  secondaries  like  PETN  are  suggestive  of  many 
other  compounds  that  exist  between  the  main  body  of  primaries  and  the 
main  body  of  secondaries.  However,  these  have  received  relatively 
little  attention  because  there  is  no  obvious  practical  use  for  them. 
Summarizing  from  a physico-chemical  point  of  view,  primaries  and  secon- 
daries differ  In  degree  rather  than  In  kind. 

The  most  important  properties  of  a primary  explosive  are  the  ability 
to  remain  ignited  and  the  ability  to  undergo  very  fast  deflagration  to 
detonation  transition  (DDT).  These  ignition  properties  are  a primary's 
strength  and  its  weakness.  It  is  a strength  because  the  function  of  a 
primary  is  to  detonate  on  being  l-inlted:  all  other  things  being  equal, 

the  faster  the  DDT,  the  better  the  primary.  It  is  a weakness  because 
accidental  initiation  of  deflagration  will  also  result  in  detonation. 
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Because  DDT  seemed  to  be  a distinguishing  characteristic  that  could 
be  measured,  an  experimental  program  was  designed  to  test  whether  DDT 
measurements  could  be  used  as  a definitive  screening  technique.  This 
research  program  was  based  on  the  fo7  -.owing  assumptions: 

(1)  An  important  measurable  distinction  between  a primary  and  a 
secondary  explosive  is  that  a primary  explosive  undergoes 
very  much  faster  deflagration  to  detonation  transition  (DDT). 

(2)  The  difference  between  primaries  and  secondaries  is  one  of 
degree  rather  than  kind;  therefore,  experience  obtained  in 
measuring  DDT  in  secondaries  can  be  applied  to  primaries. 

(3)  It  would  be  easier  to  develop  a quantitative  model  for  DDT 
in  cast  explosives  than  in  porous  explosives. 

DDT  studies  of  two  potential  isomeric  cast  primaries,  l-methyl-5- 
nltrotetrazole  (l-M-'IT)  and  2-methyl-5-nitrotetrazole  (2-MNT),  were 
undertaken  as  the  first  step  towards  achieving  this  objective.  Isomeric 
cast  explosives  are  more  convenient  for  isolating  the  effects  of  chemical 
structure  in  ( rimaries  because  isomers  have  different  molecular  structures 
but  essentially  the  same  physical  properties;  because  inhomogeneities 
leading  to  increased  sensitivity  are  minimized  by  the  casting  process; 
and  because  the  mechanism  of  DDT  is  more  fully  understood  in  cast 
explosives  than  in  porous  explosives. 

Our  experiments  showed  chat  DDT  occurs  more  readily  in  cast  2-MNT 
than  in  cast  1-MNT,  clearly  demonstrating  the  influence  of  chemical 
structure  on  primary  explosive  behavior.  The  threshold  voltages  for 
ignition  by  a hot  bridgewire  were  found  to  be  110  volts  for  cast  2-MNT 
and  375  volts  for  cast  1-MNT.  Critical  times  and  distances  for  DDT  in 
cast  2-MNT  of  8 fxsec  and  2-3  mm,  respectively,  were  determined  by 
performing  experiments  on  charges  of  different  lengths.  The  same  critical 
parameters  for  cast  1-MNT  were  not  determined,  as  will  be  explained  later, 
because  the  threshold  voltage  for  ignition  was  too  high. 
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Macek's  model  for  DDT  in  cast  explosives  was  modified  in  the 
theoretical  study  to  provide  a more  realistic  treatment  of  the  accelerating 
flame.  Calculations  were  made  with  the  deflagration  treated  as  a reactive 
discontinuity  whose  burnt  state  satisfied  either  the  Chapman-Jouguet  (CJ) 
condition  or  the  condition  of  zero  particle  velocity  to  determine  upper 
and  lower  limits  for  the  deflagration  velocity  as  functions  of  pressure 
in  the  unburned  explosive.  Times  and  distances  for  DDT  in  cast  pentollte 
and  crystalline  lead  azide  were  predicted  and  agreed  well  with  the  meager 
experimental  data  available. 


II  STRUCTURAL  CONSIDERATIONS 


A.  Structure-Property  Characteristics 

The  objective  of  the  first  part  of  the  general  program^  was  to 
survey  the  published  literature  to  find  correlatable  characteristics 
of  primary  compounds.  Complementary  to  this  effort,  we  interviewed  a 
number  of  research  workers  familiar  with  primary  explosives  to  gather  the 
latest  unpublished  information  on  compounds  other  than  azides  under 
development.  Their  views  were  also  solicited  on  which  aspects  of  research 
were  most  needed  to  provide  new  candidate  compounds  for  useful  detonants. 

Emphasis  in  the  survey  and  discussions  was  on  compounds  that  were 
generally  considered  "primaries"  by  their  sensitivity  characteristics 
or  usage.  Table  1 lists  thirteen  general  classes  or  families  of  compounds 
and  their  principal  energetic  groupings.  The  literature  was  searched  for 
structure  characteristics  that  may  be  a key  to  their  "primariness;" 
thermochemical  properties  such  as  heat  of  formation;  reaction  kinetics; 
detonation  properties,  including  detonation  velocity,  pressure,  ease 
of  initiation  by  hot  wire,  deflagration  to  detonation  information;  and 
physical  properties  of  density,  stability,  sensitivity,  crystal  forms, 
and  related  data. 

The  survey  did  not  reveal  a consistently  measured  property  that 
would  allow  a correlation  in  a broad  sense  for  predictive  purposes, 
except  for  the  commonly  experienced  sensitivity  characteristics.  The 


Not  all  compounds  considered  as  primary  explosives  are  useful  as 
initiators  or  detonants;  the  latter  are  those  compound."?  that  have  passed 
a number  of  practical  tests  and  development  in  the  practical  configu- 
ration of  a detonator  or  fuze  train.  An  expression  sometimes  heard  is, 
"A  primary  the  compound  may  be,  but  a detonant  it  is  not." 


6 


Table  1 


Classes,  Energy  Group 


Acetonl trlles 


-C  - CN 
Acetylldes 


-C  s CH 


Azides 

-Na 

Fulminates 

-NCO 

Furoxanes 


O 

N, 


/ 


/ 


N1  trocyanatnides 

-N(NOj)CN 

Nitroaromatics 


(NOa 


CLASSES  OF  PRlKARy  EXPLOSIVES 


Salts 


(OH) 


K-dlnitroacetonile 


Ag,  Ba 


Pb  Styphnate 


K,  Ba,  Trinltrophloro- 
gluclna te 

K,  Ba,  Dinitrocatechol 


Covalent 


Ag 


Pb,  Ag,  Cu,  Th 


Hg,  Ag  , other.'- 


K dinltrobenzofuToxane  Benzotrifuroxan 


Trlnitroacetonitrlle 

Dlcyanoacotylene 

Trini trotrlazidobenzene 
Cyanuric  triazide 


Ni trocyanamlde  *H  jO 

Stypbnic  Acid 

Trini trophloroglucinol 

Dinl trocatechol 

Dlazodinitrophenol , R=Nj 

Bi s(d if luor amino )hexa- 
nitrodiphenyl  (FDIPAM), 
R=KF, 


Table  1 (Concluded) 


CLASSES  OF  PRIMARV  EXPLOSIVES 


Classea,  Energy  Group 


Salts 


Tatrazoles 


N N 


\ 

0 


,C  - R 


R = H(Ag,  Cu,  Ni,  Pb) 

R = NHs,(Nl,  Co,  Pb) 

R = N3<Ag,  Na,  K,  Cd^ 

R = NN-tetrazole(Hg,  Tl,  Na) 
R = CKAg.  Cu) 


R = NH-NH-tetrazolef Na,  Ag, 
Tl,  Cu,  Hg,  Ni,  Co,  Pb) 

R = CH,(Ag,  Cu,  Pb) 

R = Phenyl (Hg,  Cu,  Ph^ 

a - NHNOj(Ag,  Tl,  K,  Pb) 

R = NO,(Hg,  Ag,  Tl,  Na,  K, 
Cu,  Hg(ic),  NI,  Co,  Pb) 

R « NHNHj(Ag,  Tl,  K,  Cu,  NI, 
Co,  Pb) 

Perchlorate  CP 


X 


Nil, 

X = N N>. 

II  C - CN 

N N ' 


(CIO^), 


Miscellaneous 


Nitrogen  triiodide 
Dini trogentetrasulfide 
Ag  - methylene  dinitramine 


8 


1.  i ■-  ' • ‘ 


Covalent 

Tetrazols  (unstable) 
Tetrazene 

1- Methyl-5-ni tro tetrazole 

2- Methyl“5-ni trotetrazole 
Azidotetra zole 


latter  will  for  some  time  probably  broadly  define  primary  candidates 
as  being  compounds  easily  Ignited  by  low  energy  inputs  from  impact, 
spark,  hot  wire,  and  shock.  Most  of  the  data  that  were  even  moderately 
useful  for  correlation  were  usually  within  families  of  compounds. 

There  is  a priori  no  basis  at  present  for  assuming  primariness  from 
structure  examination.  However,  one  may  Identify  energetic  structural 
elements  in  a given  highly  energetic  family  of  compounds  that  would 
forewarn  one  to  look  for  primary  characteristics. 

For  stability  and  ease  of  handling,  the  most  used  compounds  are 
the  salts  of  highly  sensitive  covalent  energetic  compounds,  such  as 
hydrazoic  acid  and  the  tetrazoles,  and  the  salts  of  some  less  sensitive 
secondary  explosive  compounds,  such  as  styphnic  acid.  The  inclusion  of 
azo  linkages  in  a structure  seems  to  impart  primary  sensitivity,  especially 
in  the  form  of  the  azide  group,  -N3,  and  in  heterocyclic  high  nitrogen 
compounds  as  represented  by  the  triazole  and  tetrazole  groups.  Similar 
unsaturate  nitrogen  linkages  are  i*epresented  by  the  fulminate  NCO  and 
by  the  hetero  furoxane  in  the  "free"  form,  such  as  benzotrifuroxane  (BTF) 
or  salt  form,  e.g.,  potassium  dinitrobenzofuroxan.  Other  primary 
candidates  have  come  from  nitrogen  unsaturate  compounds  having  combina- 
tions of  the  above  groups,  as  represented  by  the  nitro  furoxan  and  nitro 
azldo  compounds;  also  in  some  cases  the  nltroamino  group  is  combined  with 
an  unsaturate  group,  as  in  a nitrocyanamide  salt.  Some  exceptions  to 
this  classification  of  generalitiec  are  the  inorganic  sensitive  compounds 
such  as  perchlorate  salts.  Also  organic  difluoramlno  compounds  have 
sensitivity  characteristics  that  at  least  suggest  they  be  considered. 

Although  most  of  the  primaries  in  practical  use  as  detonants  are 
salts,  an  increasing  number  of  covalent,  or  "free  acid,"  forms  of 
explosives  are  receiving  attention.  The  compound,  FDIPAM,®  was  recently 
evaluated  as  a primary.  It  was  found  to  be  comparable  to  lead  styphnate 
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In  developing  Initiation  to  detonation  and  slightly  less  sensitive  to 
Initiation  by  hot  wire  than  lead  azide.  ERDA  is  currently  evaluating 
the  compound,  BTF,  as  an  Initiator.  In  our  laboratories  tests  on  2- 
methyl-5-nltrDtetrazole  have  Indicated  that  It  has  a deflagration  to 
detonation  transition  similar  to  lead  azide;  further  study  may  be 
warranted.  Other  organic  compounds  that  have  been  considered  are  tri- 
nitrotrla^^ldobenzene,  cyanurlc  trlazlde,  and  dlazodlnltrophenol ; the 
last  has  use  as  an  Initiator  In  commercial  compositions. 

B,  Structure-Property  Correlations 

1 . Metal  Ionization  Potentials* 

A few  structure-property  correlations  have  been  suggested  for 
primary  explosives.  One  of  the  most  interesting  was  proposed  independently 
by  Gray'*  and  by  Taylor  and  Jenkins®  (GTJ).  The  GTJ  postulate  is  that, 
within  a family  of  explosive  salts  such  as  the  azides  or  tetrazoles,  the 
greater  the  Ionization  potential,  of  the  cation,  the  more  likely  the 
compound  is  a primary  explosive.  This  concept  may  be  useful  for  screening 
new  families  of  salts.  Since  H and  Ag  have  high  ionization  potentials, 
if  the  free  acid  or  silver  salt  of  a new  family  does  not  show  primise,  it 
Is  likely  that  other  salts  will  not  be  useful.  The  basis  of  this  cor- 
relation is  probably  the  variation  in  heat  of  formation  of  the  salt  with 
the  Ionization  potential  of  the  cation. 

For  example,  consider  the  azides.  Azide  thermochemistry  is  par- 
ticularly simple  because  the  main  mode  of  decomposition  is  to  the 
elements;  therefore,  the  heat  of  detonation  (an  indicator  of  sensitivity) 
is  simply  the  heat  of  formation  for  the  metal  azides, 

AgN  -♦  Ag  + 3/2  Na  + ^^(AgNg) 

and  for  hydrogen  azide, 

HNg  1/2  Hj  + 3/2  Na  + tfl^CHNg) 
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Figures  1 and  2 show  that  the  heats  of  formation  correlate 


well  with  the  Ionization  potential  of  the  metal.  The  Ionization  potential 
gives  a quantitative  measure  of  the  degree  of  Ionic  character.  It  Is 
Interesting  to  note  that  ZnNg  Is  significantly  off  the  curve  for  divalent 
cations.  Perhaps  the  heat  of  formation  of  ZnNg  should  be  remeasured  as 
the  Ionization  potential  Is  surely  well  known. 

Although  the  dTJ  postulate  Is  useful  for  screening,  the  concept 
Is  limited  because  some  exceptions  occur  and  because  It  cannot  be  used 
outside  the  family  of  compounds.  Furthermore,  the  heat  of  detonation 
Is  not  the  sole  criterion  for  usefulness  as  a primary.  TNT  has  a higher 
heat  of  detonation  than  lead  azide. 

2.  Bond  Rearrangement^ 


j 

J- 

r 


Roberts  has  postulated®  that,  for  a compound  to  be  a primary 
explosive,  the  structure  should  be  capable  of  decomposing  to  products 
with  rearrangement  of  the  electrons  (l.e.,  bonds)  rather  than  the  atoms. 
This  Is  a stimulating  Idea  and  readily  explains  the  fast  DDT  of  the 
simple-molecule  primaries  like  the  azides,  fulminates,  and  acetylldes. 


We  selected  two  compounds  for  further  study: 


J 

i 

r N = CHv 

I I /C  = 0 and 

; N = CH' 


A 


HC 

I 

HC 


= N 


= N 


> 


B 


O 


Both  isomers  formally  dissociate  with  only  electronic  and  not  atomic 
rearrangement  to  2HCN  + CO.  (There  are  no  reports  In  the  literature  of 
synthesis  of  A or  B and  no  new  routes  to  A and  B that  seem  reasonable.) 
Benson  and  Shaw  have  estimated^  the  heat  of  formation  of  A to  be  43  ± 

5 kcal/mole  and  of  B to  be  1 ± 5 kcal/mole. 


I 
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SA-2842-1 


FIGURE  1 CORRELATION  OF  HEATS  OF  FORMATION  IN  THE  SOLID  STATE 

OF  MONOVALENT  AZIDES  WITH  FIRST  IONIZATION  POTENTIAL  OF 
THE  CATIONS 
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FIGURE  2 CORRELATION  OF  HEATS  OF  FORMATION  IN  THE  SOLID  STATE  OF 
DIVALENT  AZIDES  WITH  SECOND  IONIZATION  POTENTIAL  OF 
THE  CATIONS 
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The  minimum  rearrangement  products  for  both  A and  B are  2HCN 
CO.  The  sum  of  the  heat  of  formation  of  the  products  Is  (2  x 32,3)  - 
26.4  a 38.2  kcal/mole.  Therefore^  decomposition  of  A to  2HCN  and  CO 
would  liberate  5 5 kcal/mole,  and  decomposition  of  B to  the  same 

products  would  adsorb  37  ± 5 kcal/mole  from  the  surroundings.  Thus  A 
Is  a marginal  explosive,  and  B Is  unlikely  to  explode  by  electron  rear- 
rangement mechanism. 

Another  argument  against  the  concept  arises  from  examining  the 
gas  generator,  (-cnaN-NsK:H-)jj,  prepared  by  Reed,®  Benson  and  Shaw 
have  estimated  the  heat  of  formation  of  the  repeating  unit  to  be  68  5 

kcal/mole.  If  the  products  ai*e  2HCN,  (2  x 32.3)  - 64.6  kcal/mole,  the 
decomposition  yields  3.4  kcal/mole,  and  If  the  products  are  CHaCH  (54.2 
kcal/mole)  -t-  N,  (0  kcal/mole.  Reed  found,  consistent  with  these  estimates, 
that  the  compounds  Indeed  did  not  detonate  and  liberated  very  little  heat 
when  they  decomposed . 

Finally,  a very  persuasive  argument  against  the  minimum 
rearrangement  concept  Is  the  structure  of  styphnlc  acid,  a primary 
explosive.  OH 

NOj 

OH 

It  Is  hard  to  see  how  that  molecule  can  decompose  without  considerable 
atomic  rearrangement, 

3.  Substltuency  Effects  In  Tetrazoles 

The  fulminates  and  tetrazoles  are  the  only  classes  of  compounds 
other  than  azides  that  have  received  systematic  Investigation  of  proper- 
ties, the  latter  principally  at  the  Explosives  Research  and  Development 
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Establishment  In  Waltham  Abbey.  As  listed  in  Table  1,  a large  number  of 
tetrazole  salts  and  covalent  structures  have  been  synthesized  and  charac- 
terized for  usefulness  as  detonants.  One»  mercury  3-nitrotetrazole,  has 
also  been  investigated  by  the  U.S.  Navy  and  is  now  proposed  as  a replace- 
ment for  load  azide.’’ 

Sufficient  data  on  tetrazoles  have  been  accumulated  to  allow 
research  workers  to  make  a few  structure/property  correlations.  Bates 
and  Jenkins^’’  have  suggested  that  the  ranking  in  explosive  behavior  of 
twelve  3-substltuted  nitrotetrazoles  could  be  related  to  the  substituent's 
electron  withdrawing  power:  the  more  electron  withdrawing,  the  more 

explosive  the  compound.  The  ranking  observed  for  the  substituent  in 
the  5 position  ranged  from  mild  ignitions  for  methyl  5-nltrotetrazole 
to  instability  in  the  following  order:  CHg  =»  CgH^  (mild  ignitions), 

< NHj  < H < NHNOj  (compounds  explode),  < bis( tetrazole)  < 5,5 '-azo- 

dl tetrazole  (do  not  detonate  RDX),  < Cl  < NOg  (very  powerful  explosives) 

< Ng  (very  sensitive)  < -Nj  (unstable). 

The  salts  of  5-nltrotetrazole  could  also  be  ranked  according 
to  their  sensitiveness.  The  salts  of  silver,  mercury,  and  lead  gave 
more  sensitive  compositions  and  greater  initiating  power  than  the  alkali 
metal  salts.  The  same  salts  of  5-chlorotetrazole  were  too  sensitive 
and  corrosive,  and  those  of  5-azidotetrazole  were  too  hazardous.'^ 

Haskins^*  has  compared  molecular  orbital  calculations  of  some 
azide  and  tetrazoles  with  explosive  behavior,  dipole  moments,  and 
binding  energies.  He  used  an  iterative  extended  Huckel  molecular 
orbital  approach  to  calculate  dipole  moments  of  l-(H)tetrazole,  2-(H)tetra- 
zole,  and  5-aminotetrazole  and  found  good  agreement  with  experimental 
measurements.  For  example,  1-(H) tetrazole  was  calculated  as  5.269^ 
and  measured  as  5.11|j.  Binding  energies  calculated  by  the  similar 
extended  extended  Huckel  molecular  orbital  (EHMO)  method  for  the  1-(H) 
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and  2-(H)  compounds  were  also  In  good  agreement  with  experimental 
measurements. 

Haskins  also  compared  tae  amplrloai  ranking  of  exploslv* 
behavior  proposed  bj'  Jenkins  with  one  devl  lod  according  to  total  charge 
In  the  substituent  group.  Thus,  EHMO  results  gave  a quantitative 
measure  of  the  electron  withdrawing  power  of  the  substituents  and  gave 
an  order:  CH,  < H < NH,  < NHNO,  < N,  < NO,.  Haskins  suggested  that  the 
EHMO  calculations,  although  not  In  agreement  with  the  qualitative 
ordering,  were  enough  In  agreement  to  Indicate  a relationship  between 
explosive  behavior  and  electron  withdrawal  power  from  the  ring. 

4.  Physical  Variables 

The  basic  problem  with  structure-property  correlations  in 
primary  explosives  Is  the  scarcity  of  quantitative  data  for  a very  large 
number  of  variables.  From  the  structure  standpoint  alone,  the  number  of 
candidate  primaries  within  the  families  listed  In  Table  1 Is  very  large, 
although  they  come  within  the  classification  generalities  discussed 
earlier.  Overlaying  all  chemical  properties  are  the  physical  variables, 
which  may  determine  whether  a primary  becomes  a detonant.  Physical 
variables  of  paartlcle  size  and  crystal  form  produce  a wide  variation 
In  applicability  to  a detonator  configuration;  lead  azide  Is  a prominent 
example  of  physical  variable  Influence.  Even  an  explosive  such  as 
hexanltrostllbene  (HNS),  generally  considered  as  secondary,  can  have 
primary  sensitivity  characteristics  in  superfine  particle  size,  HNS  I, 
in  contrast  to  HNS  form  II  which  is  relatively  insensitive.^®  It  appears 
that  none  of  the  sensitivity  data  can  be  used  as  a qualifying  teat  for 
a primary. 

However,  we  have  proposed  that  one  generally  unmeasured  property 
of  primaries,  deflagration  to  de'onation  transition  (DDT),  can  be  used 
for  screening  and  classifying  primary  compounds,  The  property  of  fast 
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DDT  (characteristic  of  lead  azide),  when  combined  with  the  observed 
characteristic  that  a primary  will  not  extinguish  when  once  ignited,  may 
be  a criterion  for  distinguishing  primaries  from  secondaries,  as  a 
function  of  structure  and  Independent  of  sensitivity  properties.'*’'  The 
distinction  at  present  Is  therefore  one  of  degree;  defining  a primary  In 
absolute  terms  1s  not  possible  with  currently  available  Information, 


Current  practices  contraindicate  the  use  of  DDT  measurements  as  a neces- 
sary criterion;  e.g.,  some  primaries  have  long  DDT  lengths,  such  as  lead 
styphnate. 
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Ill  DEFLAGRATION  TO  DETONATION  IN  CAST  EXPLOSIVES 
A.  Background 

For  DDT  measurement  to  be  a practicable  screening  procedure,  a model 
system  Is  needed  to  establish  the  basic  approach.  The  objectives  of  the 
work  reported  here  were  to  determine  If  such  a measurement  would  permit 
one  to  relate  the  prlmarlness  of  an  experimental  chemical  structure,  to 
a DDT  result,  In  comparison  with  DDT  test  results  of  practical  compounds 
such  as  the  azides.  Chemical  structure  relationships,  primarily  the 
effect  of  substituents,  can  be  obtained  by  measuring  a number  of  analog 
molecules  of  a class  (as  in  the  tetrazoles)  or  by  measuring  the  effect 
of  structural  changes  in  isomers.  The  latter  allows  deductions  to  be 
made  on  compounds  having  the  same  empirical  formula  but  differing  in 
the  position  of  a substituent. 

A second  requirement  for  this  study  was  that  the  model  compounds  be 
castable  since  much  DDT  state  of  the  art  has  been  developed  with  cast 
secondaries.  A survey  of  known  "low  melting"  primaries  revealed  that 
both  the  1-  and  2-methyl  derivatives  of  5-nitrotetrazole  fulfilled  these 
requirements. 

The  research  plan  was  based  on  the  following  assumptions: 

(1)  An  important  distinction  between  a primary  and  a secondary 
explosive  is  that  a primary  explosive  undergoes  very  much 
faster  deflagration- to-detonation  transition  (DDT), 

(2)  The  difference  between  primaries  and  secondaries  is  one  of 
degree  rather  than  kind;  therefore,  experience  obtained  in 
measuring  DDT  in  secondaries  can  be  applied  to  primaries. 

(3)  It  would  be  easier  to  develop  a quantitative  model  for  DDT 
in  cast  explosives. 
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The  first  assumption  is  based  on  the  prime  functions  of  primary  and 
secondary  explosives.  Primaries  are  required  to  detonate  rapidly  when 
thermally  ingited,  whereas  detonations  in  secondaries  are  initiated  by 
shock  waves.  Possibly  the  faster  the  DDT,  the  better  the  primary  although 
this  premise  requires  a larger  data  base  than  is  now  available  to  support 
it.  The  inevitable  DDT  of  a primary  is  also  its  weakness  because  ac- 
cidental initiation  of  deflagration  will  result  in  detonation. 

The  second  assumption  is  based  on  observations  of  the  many  similarities 
between  primary  and  secondary  explosives.®*^  Both  primaries  and  secondaries 
will  undergo  DDT,  although  the  transition  takes  longer  and  requires  heavier 
confinement  for  secondaries.  The  existence  of  near  primaries,  such  as 
some  of  the  hydroxynitroaromatic  compounds,  and  sensitive  secondaries  like 
PETN  suggests  that  there  may  be  many  other  compounds  between  the  main 
body  of  primaries  and  the  main  body  of  secondaries  that  have  received 
little  attention  because  there  is  no  obvious  practical  use  for  them. 

The  third  assumption  is  based  on  the  results  of  DDT  research  on 
secondary  explosives.  A first-generation  model  that  gives  order  of 
magnitude  agreement  with  experimentally  determined  times  and  distances 
for  DDT  in  cast  secondaries  has  been  formulated  by  Macek.^‘‘  No  equivalent 
model  of  DDT  in  porous  secondaries  exists.  The  work  of  Bernecker  and  Price'® 
indicates  that  the  formation  of  a convective  flame  front  over  a large 
surface  area  and  a resulting  complex  flow  field  preceed  DDT  in  a porous 
explosive-based  propellants,  a clearer  understanding  of  the  role  of 
convective  burning  in  DDT  in  proous  explosives  may  eventually  be  developed. 

At  present,  however,  the  basic  physical  processes  involved  in  DDT  have 
been  determined  only  for  cast  explosives  in  which  convective  burning 
does  not  occur. 

Many  of  the  physical  variables  found  in  crystalline  test  compounds, 
such  as  density  and  crystal  size  and  shape,  were  eliminated  by  using  cast 
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primaries.  The  most  important  variable,  chemical  structure,  was  retained 
by  studying  two  isomers,  l-methyl-5-nitrotetrazole  (1-MNT)  and  2-methyl- 
5-nitrotetrazole  (2-MNT).  These  compounds  were  synthesized  as  described 
in  Appendix  A,  and  their  heats  of  formation,  detonation  velocities,  and 
Chapman- Jouguet  pressures  were  calculated  (see  Appendix  B) . 

The  specific  objectives  were  to  measure  the  times  to  DDT  for  cast 
1-MNT  and  2-MNT  and  to  develop  and  teat  a first-generation  model  of 
DDT  in  cast  primaries  that  relied  on  experimental  input  data.  To  meet 
these  objectives,  the  work  was  divided  into  the  following  five  tasks. 


Task  1:  Measure  the  time-for-DDT  on  cast  1-  and  2-methyl- 

5-nitrotetrazole  (1-MNT  and  2-MNT). 

Task  2:  Conduct  well-instrumented  hot-wire  initiation 

experiments  on  1-MNT  and  2-MNT  to  measure  the 
pressure-time  histories  in  the  deflagration  regions 
and  the  times  and  distances  for  detonation. 


Task  3;  Analyze  Macek's  model  for  DDT  in  cast  secondary 

explosives  and  apply  it  to  DDT  in  cast  primaries. 

Task  4;  Test  this  DDT  model  on  the  existing  DDT  data  for 
single  crystals  of  lead  azide. 

Task  5:  Use  the  experimental  pressure-time  data  obtained  in 

Task  2 to  predict  the  times  and  distances  for  DDT 
for  1-MNT  and  2-MNT  and  compare  the  predictions  with 
the  experimental  results. 


B.  DDT  Measurements  on  Cast  1-MNT  and  2-MNT 


The  time  for  DDT  was  measured  for  cast  1-MNT  and  2-MNT  by  using 
the  technique  developed  by  Leopold^®  and  modified  by  Goodale.®  The 
initiator  plug  and  charge  holder  are  shown  in  Figure  3.  The  bridgewires 
were  1-mil  Nichrome  wires  soldered  across  the  center  of  the  charge  holder. 
A plastic  monofilament  light  pipe  of  40  mils  diameter  was  placed  in  the 
phenolic  plug  with  its  upper  surface  immediately  under  the  bridgewire 
to  determine  the  time  of  initiation  of  burning  in  the  explosive  by 
observing  the  first  light  emission  near  the  bridgewire.  A second  light 
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pipe  was  Inserted  Into  the  top  of  the  charge  holder  and  placed  directly 
on  the  upper  surface  of  the  cast  explosive  to  detect  the  light  emitted 
by  the  detonation  wave  as  it  reached  the  top  of  the  charge. 

The  capacity-discharge  firing  circuit  diagrammed  in  Figure  4 was 
used  to  ignite  the  explosive  charges.  The  waveforms  of  the  current 
discharge  through  the  bridgewire,  as  represented  by  the  voltage  drop 
across  a constant,  noninductive  resistance  close  to  tho  bridgewlre  in 
the  circuit,  was  monitored  in  each  test  on  an  oscilloscope.  Reproducible 
current  pulses  free  from  anomalous  irregularities  required  the  use  of  a 
vacuum  switch  of  high  current  capacity.  The  resistance  of  each  bridgewlre 
was  measured  before  each  test  to  ensure  that  the  bridgewlre  connections 
remained  Intact  and  that  the  resistance  remained  in  the  normal  range 
of  4.5  to  7.5  ohms. 

The  electronic  recording  system  was  slightly  different  from  that 
used  previously,®  As  in  the  previous  study, a dual-beam  oscilloscope 
was  used  to  record  the  capacitor-discharge  pulse  and  the  light  emitted 
by  the  ignition  of  the  explosive.  A single-beam  oscilloscope  with  a 
very  sensitive  triggering  mechanism  was  used  to  record  the  light  output 
observed  by  the  lower  light  pipe.  A pulse  generator  was  added  to  the 
circuit  to  p'"oduce  a time  reference  peak  at  0.6  ^sec  after  the  triggering 
of^the  lower  photomultiplier  by  the  light  emitted  by  the  ignition  of 
the  charge.  To  facilitate  observation  of  the  entire  trace,  the  triggering 
pulse  was  delayed  for  1.57  ^sec  and  then  sent  to  the  lower  photomultiplier. 
Another  single-beam  oscilloscope,  also  triggered  by  the  ignition  of  the 
explosive,  was  used  to  record  the  light  output  from  the  upper  light  pipe. 

Accurately  weighed  5-,  10- , 20-,  and  40-mg  charges  of  1-MNT  and 
2-MNT  were  transferred  into  charge  holders  and  pressed  at  5000  psi. 

These  pressed  charges  were  then  melted  to  obtain  cast,  fully  dense 
charges.  The  initial  castings  of  2-MNT  contained  numerous  voids  and 
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C - CaDacitor  - t.O  microfarad 
S - Vacuum  Switch 

I.  - Circuit  Inductance  - 1 6 microhenries 
- Circuit  Ftesistence  - 0.31  ohm 
Rj  - Coaxial  Current  Shunt  - 0.006  ohm 
Rg  - Bridgewire  - 4.5  to  7.5  otiins 

SA-4770-1 

FIGURE  4 FIRING  CIRCUIT 


rough  surfaces  because  they  were  allowed  to  solidify  too  rapidly.  A 
casting  procedure  was  developed  in  which  a slow  cooling  rate  and  contin- 
uous observation  through  a microscope  were  used  to  ensure  reliable  casting. 
The  resulting  cast  2-MNT  charges  had  extremely  flat  upper  surfaces  and 
an  average  density  of  1,672  g/cm®.  Considering  the  small  weights  (5-40  mg) 
and  measured  charge  lengths  (0.5-5  mm),  this  average  density  was  in  good 
agreement  with  the  literature  value^’  of  1.64  g/cm®. 

The  same  casting  procedure  was  then  used  to  prepare  the  cast  1-MNT 
charges.  Because  of  its  low  melting  point  (58^C  compared  with  86°C  for 
2-MNT),  1-MNT  cooled  relatively  slowly  and  exhibited  a large  degree  of 
supercooling  before  solidifying  Into  a glass.  By  using  a faster  cooling 
rate,  fully  dense,  crystalline  charges  of  1-MNT  were  finally  produced  in 
the  forms  of  microcrystals  and  relatively  coarse  crystals.  All  these 
forms  of  cast  1-MNT  (glassy,  microcrystalline,  and  coarse  crystalline) 
were  tested.  The  average  density  of  the  cast  1-MNT  charges  (1.775  g/cm^) 
closely  agreed  with  the  value  of  1,76  g/cm®  predicted  by  a group  additi- 
vity approach  previously  described.* 

Since  both  cast  1-MNT  and  2-MNT  were  found  to  be  somewhat  transparent 
to  light,  the  tops  of  all  charges  were  blackened  with  a thin  film  of 
India  ink  to  prevent  premature  triggering  of  the  upper  photomultiplier 
by  light  transmitted  from  inside  the  reacting  charge  before  the  detonation 
wave  reached  the  top  surface  of  the  charge. 

To  test  the  apparatus,  we  performed  a series  of  time-for-DDT  measure- 
ments on  pressed  lead  azide  charges  prepared  four  years  ago,*  The  mean 
time  to  DDT  for  these  charges  was  0,9usec  as  compared  with  0.7  fisec  ob- 
tained previously.  However,  the  range  of  deviations  was  similar  to  the 
range  of  deviations  observed  for  freshly  prepared  lead  azide.  The  slightly 
longer  times  to  DDT  may  have  been  due  to  aging  effects  in  che  lead  azide. 
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Two  series  of  cast,  20-mg  charges  were  fired  at  various  capacitor 
charging  voltages  to  determine  the  threshold  voltages  for  initiation 
of  1-MNT  and  2-MNT,  For  cast  2-MNT,  the  threshold  voltage  for  20-mg 
charges  was  approximately  100  volts,  while  discharges  of  110  volts 
reliably  detonated  every  2-MNT  charge  fired.  Therefore,  capacitor 
charging  voltages  of  110  volts  were  used  to  make  the  quantitative 
measurements  of  time  to  DDT  as  a function  of  charge  length  for  cast 
2-MNT.  Although  it  was  expected  that  cast  1-MNT  would  have  a similar 
threshold  voltage,  cast  1-MNT  charges  failed  to  explode  at  capacitor 
charging  voltages  of  less  than  375  volts,  and  only  one-fourth  of  the 
charges  exploded  at  that  voltage.  In  most  tests,  sufficient  pressure 
was  generated  to  blow  the  charge  holder  off  the  plug.  Usually  unreacted 
1-MNT  remained  in  the  charge  holder,  and  no  audible  explosive  effects  were 
produced.  At  capacitor  charging  voltages  above  200  volts,  the  light 
emitted  by  the  glowing  bridgewire  was  sufficient  to  trigger  the  sensitive 
lower  photomultiplier  before  the  explosive  actually  Ignited.  Thus  the 
[ measured  times  to  DDT  for  these  high  voltage  tests  were  upper  limits 

[ because  they  Included  the  ignition  delay  time  of  cast  1-MNT. 

I To  avoid  this  premature  triggering  at  high  voltages,  a modified  bridge- 

! wire  apparatus  was  fabricated  in  which  the  bridgewire  was  moved  off  center 

I 

I so  that  its  glow  would  not  be  as  readily  sensed  by  the  lower  light  pipe. 

I This  modified  apparatus  worked  well  at  voltages  below  approximately  350 

|.  volts,  but  bridpewire  glow  remained  a problem  at  the  voltage  required 

I f ' r cast  1-MNT  initiation.  Therefore,  quantitative  time  to  DDT  versus 

! cnarge  length  measurements  could  not  be  obtained  for  cast  1-MNT. 

However,  the  large  difference  in  threshold  capacitor  charging 
voltages  for  two  cast  isomers,  1-MNT  and  2-MNT,  is  an  important  result 
; bci-ause  it  indicates  that  chemical  structure  can  play  a major  role  in 

the  performance  of  potential  primaries.  Because  the  threshold  voltage  of 
[ cast  2-MNT  is  approximately  the  same  as  that  of  pressed  lead  azide, ^ 


cast  2-MNT,  based  solely  on  this  test,  appears  to  be  a promising  candidate 
primary  explosive,  while  cast  1-MNT  appears  not  to  be  a primary.  At 
present  the  only  tentative  explanation  of  this  behavior  is  that  in  1-MNT 
the  close  proximity  of  the  hydrogen  atoms  in  the  methyl  groups  to  the 
oxygen  atoms  in  the  nitro  group  leads  to  sufficient  hydrogen  bonding  and 
decreased  sensitivity. 

Quantitative  data  on  time  for  DDT  as  a function  of  charge  length 
were  obtained  for  cast  2-MNT  using  a capacitor  charging  voltage  of  110 
volts.  These  cast  charges  varied  in  weight  from  4 to  47  mg  and  in  length 
from  0.4  to  5 mm.  Charges  weighing  less  than  5 mg  failed  to  detonate, 
while  approximately  two- thirds  of  the  10-mg  charges  detonated.  Charges 
exceeding  20  mg  definitely  detonated  and  damaged  the  aluminum  charge  holder 
to  roughly  the  same  extent  as  an  equivalent  mass  of  pressed  lead  azide. 

The  experimental  results  for  the  fifteen  cast  2-MNT  charges  that  produced 
good  oscilloscope  records  of  three  of  these  shots  (shots  2M-45,  2M-34, 
and  2M-19  in  Table  2)  are  reproduced  as  Figures  5,  6,  and  7,  respectively. 
Figures  5a,  6a,  and  7a  are  the  dual-beam  oscilloscope  records  of  the 
bridgewire  current  waveform  and  the  light  produced  by  the  Ignition  of 
cast  2-MNT  near  the  bridgewire.  Practically  all  the  energy  from  the 
capacitor  discharge  was  deposited  in  the  bridgewire  in  4-6  Usee,  and 
ignition  occurred  11-16  fisec  after  the  discharge  began.  The  average  time 
from  the  beginning  of  the  energy  pulse  to  the  onset  of  burning  for  the  cast 
2-MNT  charges  listed  in  Table  2 was  13.7  usee.  The  corresponding  time 
to  ignition  for  pressed  lead  azide  at  105.9  volts  was  about  50  u^ec.^ 
Figures  5b,  6b,  and  7b  are  the  single-beam  oscilloscope  records  of  the 
light  output  detected  by  the  lower  photomultiplier  after  it  had  been 
triggered  by  the  first  light  of  ignition  near  the  bridgewire.  As  pre- 
viously mentioned,  a reference  pulse  with  a 0.6  ^sec  time  delay  and  a 
time  delay  circuit  of  1.57  yisec  for  the  txace  allowed  observation  of 
the  entire  lower  photomultiplier  output  with  a known  time  reference. 
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Table  2 

SUMMARY  OF  TIME  TO  IGNITION  AND  TIME  FOR  DDT 
EXPERIMENTS  ON  CAST  2-MNT 


Charge 

Number 

Ci.pacltor 

Charging 

Voltage 

(volts) 

Charge 

Length 

(mm) 

Time  to 
Ignition 
of  Charge 
(usee) 

Time  from  Ignition 
to  Detection  of  Light 
at  Top  of  Charge 
(usee) 

Average 
Burning  Rate 
(mm4l  sec) 

2M-16 

110 

0.488 

12.8 

4.4 

0.111 

2M-27 

110 

0.632 

>20.0 

2.6 

0.243 

2M-26 

110 

0.904 

16.0 

2.4 

0.377 

2M-24 

110 

1.214 

19.0 

4.0 

0.3038 

2M-10 

110 

1.621 

12.4 

2.0 

0.8108 

2M-14 

110 

1.717 

12.5 

2.0 

0.859 

2M-33 

150 

1.796 

>10.0 

3.6 

0.499 

2M-45 

110 

1.923 

11.2 

2.7 

0.712 

2M-3S 

110 

1.991 

13.1 

3.2 

0.622 

2M-34 

100 

2.042 

13.5 

6.3 

0.324 

2M-46 

110 

2.223 

12.2 

1.9 

1.11 

2M-37 

110 

4.270 

15.8 

8.8 

0,485 

2M-19 

110 

4.425 

11.8 

8.6 

0.514 

2M-38** 

13  5 

4.493 

7.5 

0.599 

2M-20 

110 

4.877 

14.3 

9.0 

0.542 

°Aftor 

record 

failing  to  ignite 
of  the  capacitor 

at  110  volts,  Charge  2M- 
dlscharge  current  pulse 

•38  was 
or  the 

ignited  at  13 > volts.  No 
time  to  Ignition  was  ob- 
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a.  BRIDQEWIRE  CURRENT  WAVEFORM 


HorlionttI  Scalt:  2 (iMc/cm 
Vartical  Scala:  0.2  volti/om 

PHOTOMULTIPLIER  MONITORING 
LIGHT  NEAR  BRIOGEWIRE 

Horiiontal  Scala:  2 tmcfom 
Vartical  Scala:  0.2  volta/om 


b.  LOWER  PHOTOMULTIPLIER  TRIGGERED 
BY  FIRST  LIGHT  NEAR  BRIOGEWIRE 

Horiiontal  Scala:  0.5  /nMc/cm 
Vartical  Scala:  0.2  voltt/cm 

REFERENCE  PEAK  AT  0.6  piac 

TIME  DELAY  OF  1.67  mmc 


c.  UPPER  PHOTOMULTIPLIER  TRIGGERED 
BY  FIRST  LIGHT  NEAR  BRIOGEWIRE 

Horizontal  Scala:  1 M>ac/cm 

Vartical  Scala:  0.2  voltt/cm 

REFERENCE  PEAK  AT  0.6  nmc 
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FIGURE  5 OSCILLOSCOPE  RECORDS  FOR  SHOT  2M-45 


a.  BRiOQEWmE  CURRENT  WAVEFORM 

Horitontal  Scalt:  2 Miac/cm 
VtrtictI  Scalt:  0.2  volts/ctn 

PHOTOMULTIPLIER  MONITORING 
LIGHT  NEAR  BRIOGEWIRE 

Horitontal  Scalt:  2 Mwc/cm 
Vtrtical  Scalt:  0.2  voitt/cm 


b.  LOWER  PHOTOMULTIPLIER  TRIGGERED 
BY  FIRST  LIGHT  NEAR  BRIDQEWIRE 

Horiiontal  Scalt:  O.S  liUtlcm 
Vtrtical  Sralt:  0.2  voItt/cm 

REFERENCE  PEAK  AT  0.6  mMC 
TIME  DELAY  OF  1.67  mmc 


C.  UPPER  PHOTOMULTIPLIER  TRIGGERED 
BY  FIRST  LIGHT  NEAR  BRIDQEWIRE 

Horizontal  Scalt:  2 fittc/cm 
Vtrtical  Scalt:  0.2  volt*/cm 

REFERENCE  PEAK  AT  0.6  (i»c 
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FIGURE  6 OSCILLOSCOPE  RECORDS  FOR  SHOT  2M-34 
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a.  BRIDQEWIR6  CURRENT  WAVEFORM 
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Horlionul  Sealt:  2 ttttelvn 
Vtrtinl  Sctit:  0.2  voitt/om 

PHOTOMULTIPLIER  MONiTORINQ 
LIGHT  NEAR  BRIDGEWIRE 

HorlionttI  Scilc  2 imvfem 
VtrtlcN  Scilti  0.2  voltt/em 


b.  LOWER  PHOTOMULTIPLIER  TRIGGERED 
BY  FIRST  LIGHT  NEAR  BRIDGEWIRE 

Horizontal  Sctia:  0,B  Mitc/cm 
Vartleal  Scala;  0.2  volta/em 

REFERENCE  PEAK  AT  0.6  »iiac 
TIME  DELAY  OF  1.67  iMC 


c.  UPPER  PHOTOMULTIPLIER  TRIGGERED 
BY  FIRST  LIGHT  NEAR  BRIDGEWIRE 

Horizontal  Scala:  1 iMtc/cm 
Vartleal  Scala;  0.2  volls/cm 

REFERENCE  PEAK  AT  0.6  Mtac 


SA-4770-4 


I FIGURE  7 OSCILLOSCOPE  RECORDS  FOR  SHOT  2M-19 

I 
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The  duration  of  the  light  production  at  the  bottom  of  the  charge  can  be 
determined  from  records  like  Figures  5b,  6b,  and  7b, 

Figures  5c,  6c,  and  7c  represent  the  single-beam  oscilloscope  records 
of  the  light  output  detected  by  the  upper  photomultiplier  after  it  had 
been  triggered  by  the  Ignition  of  the  charge.  The  0.6-|Ltsec  reference 
peak  also  appeared  on  these  records,  which  measured  the  time  between 
Ignition  and  the  appearance  of  light  at  the  top  of  the  charge.  This 
light  was  assumed  to  be  caused  by  a detonation  wave  produced  inside  the 
charge  that  propagated  through  the  remainder  of  the  charge.  The  measured 
time  was  then  taken  as  the  time  for  DDT  for  that  charge.  The  times  for 
DDT  in  Figures  5c,  6c,  and  7c  are  2.7,  6.3,  and  8,3  |Lisec,  respectively. 
Because  shots  2M-45  (Figure  5)  and  2M-34  (Figure  6)  were  nearly  the  same 
charge  length,  the  large  difference  in  the  times  for  DDT  for  these  two 
shots  indicated  that  the  duration  of  burning  before  the  onset  of  detonation 
in  cast  2-MNT  charges  listed  in  Table  2 were  found  to  be  4.60  and  2,68  Usee, 
respectively.  These  values  are  compared  with  previously  determined  values® 
for  three  other  primaries  in  Table  3.  All  the  primaries  tested  in  this 
apparatus  exhibited  significant  standard  deviations  in  time  for  DDT.  The 
absence  of  confinement  on  the  upper  surface  of  the  explosive  may  result 
in  upward  movement  of  the  charge  during  deflagration.  This  motion  could 
contribute  to  the  variation  in  burning  times.  Cast  2-MNT  exhibits  a much 
larger  standard  deviation  in  the  times  for  DDT  than  the  other  three 
primaries,  However,  its  mean  time  for  DDT  is  less  than  that  of  lead 
styphnate  and  FDIPAM  and,  therefore,  cast  2-MNT  qualifies  as  a good 
primary  based  on  relative  DDT  times. 

If  the  course  of  DDT  were  exactly  repeatable  in  successive  experimental 
charges  of  varying  lengths,  a plot  of  charge  length  as  abscissa  against 
time  between  ignition  at  one  end  of  the  charge,  and  appearance  of  light 
at  the  other  would  appear  as  a single  well  defined  curve.  This  curve 
would  rise  from  the  origin  with  an  initially  steep  slope  representing  the 
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Table  3 


TIME  FOR  DDT  FOR  CAST  2-MNT  AND  OTHER  PRIMARY  EXPLOSIVES 


Number  of  charges  fired 

Capacitor  voltages  (volts) 

Average  charge  length  (mm) 

Mean  times  from  Ignition  to 
detection  of  light  at  the  top 
of  the  charge  (usee) 

Standard  deviations  of  times 
from  ignition  to  detection  of 
light  at  the  top  of  the  charge 
(jUsec) 


Lead 

Azide 

Lead 

Styphnate 

FDIPAM 

Cast 

2-MNT 

20 

29 

27 

15 

105.9 

67.8 

131.2 

110.0 

1.43 

1.94 

2.74 

2.31 

0.65 

11.20 

9.81 

4.60 

0.115 

1.110 

1,660 

2.68 
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initial  slow  propagation  of  the  deflagration  wave.  The  slope  would 
decrease  with  time  until  it  reached  a constant  value  corresponding  to  the 
detonation  velocity  of  the  explosive.  The  time  at  which  the  detonation 
velocity  was  reached  would  be  the  time  for  DDT  of  the  explosive. 

The  experimental  times  for  DDT  listed  in  Table  1 for  cast  2-MNT  are 
plotted  as  a function  of  charge  length  in  Figure  6.  The  times  for  DDT 
vary  greatly  for  charge  lengths  of  less  than  2 mm,  while  the  times  for 
DDT  for  4-5  mm  charges  appear  to  be  more  consistent.  The  longer  charge 
(Shot  2M-38),  which  had  a relatively  short  time  for  DDT  (7.5  |L(sec), 
failed  to  ignite  at  110  volts  and  was  fired  a second  time  at  135  volts. 

This  double  firing  probably  affected  the  performance  of  this  shot. 

Omitting  this  result,  a straight  line  with  a slope  of  3 mm/ usee  is  drawn 
in  Figure  6 through  the  results  for  the  remaining  three  longer  charges. 

The  detonation  velocity  of  cast  2-MNT  is  considerably  higher  than  3 mm/ 
jisec  (a  TIGER  code  calculation*^  resulted  in  a detonation  velocity  of 
7.8  mm/|isec),  but  the  scatter  in  the  measured  times  for  DDT  prohibits 
an  accurate  determination  of  the  detonation  velocity.  It  is  unrealistic 
to  expect  a good  value  of  detonation  velocity  from  these  experiments  on 
short  charges  because  of  the  complex  nature  of  DDT.  To  obtain  useful 
detonation  velocity  data,  the  charges  would  have  to  be  of  sufficient 
length  such  that  the  time  duration  of  the  passage  of  the  detonation  wave 
is  large  compared  with  the  time  required  for  DDT. 

In  addition  to  the  3 mm/|isec  line  in  Figure  6,  we  can  draw  another 
line  with  a much  less  steep  slope  through  the  shorter  ch.  rge  length  results. 
A line  with  a slope  of  0.1  mm/|isec  is  shown  in  Figure  6,  although  the 
large  variation  in  times  for  DDT  for  these  charges  prohibits  any  quantita- 
tive conclusions.  In  fact,  a line  with  a slope  of  0.5  mm/4sec  (the 
dashed  line  in  Figure  6)  can  be  drawn  through  all  of  the  data.  Howv’ver, 

•r  ’ 

three  observations  indicate  that  the  time-for-DDT  behavior  of  cast  2-MNT 
is  governed  by  a time-for-DDT  versus  charge-length  relationship  with  an 


abrupt  change  in  ilope  (such  as  the  solid  line  in  Figure  6).  First,  cast 
2-MNT  charges  of  appro|^\ately  1-mm  length  sometimes  failed  to  detonate, 
Indicating  that  thej  «,  r.>iical  length  for  full  development  of  detonation  is 
more  than  1 mm.  Second,  rf^arges  between  1.6  and  2 mm  did  show  a trend 
toward  large  increased  in  time  for  DDT  over  this  small  range  of  lengths. 
Third,  the  three  reliaole  sets  of  data  for  4-5  mm  charge  lengths  showed 
the  relatively  small  deviatiion  (S.8  ± 0.2  ^sec)  typical  of  a fully  devel- 
oped detonation  wave  that  bad  traveled  some  distance. 


These  three  observations  lead  to  the  conclusion  that  the  critical 
charge  length  for  DDT  is  •‘'.pproximately  2 mm  in  this  experiment.  This 
critical  distance  agre^is  qua^jj^tatively  with  two  results  for  lead  azide. 

In  a thin  flyer-plate  impact  study  of  shock- to-detonation  transition  (SDT) 
in  pressed  lead  azide,  Da, vies  et  al. found  that  the  transition  occurred 
between  1 and  2 mm  for  shock  waves  that  did  not  cause  Immediate  detonation 

at  the  explosive  surface,  ('ll>.t  did  cause  buildup  to  detonation  within  the 

■ 

explosive  sample.  Chaadhri  and  Field  concluded  ti.at  single  crystals  of 
lead  azide  less  th^-.  2 n . thick  deflagrated  rapidly  but  could  not  make  the 

..f 

transition  to  detonat|.on  that  thicker  crystals  made.  The  determination 
of  this  critical  charge  length  is  important  to  the  design  of  the  pressure 
history  experiments. 


The  main  conclusions  derived  from  the  results  of  the  time-for-DDT 
experiments  on  cast  1-MNT  and  cast  2-.MNT  are: 

(1)  Cast  2-MNT  is  a much  better  candidate  primary  than  cast 
1-MNT  because  2-MNT  charges  were  readily  initiated  by  a 
capacitor  charging  voltage  of  110  volts,  whereas  1-MNT 
charges  were  rarely  initiated  even  at  375  volts. 

(2)  Cast  2-MNT  is  a good  primary  in  its  Ignition  and  DDT 
properties  since  its  threshold  voltage  is  close  to  that 
of  pressed  lead  azide,  its  time  to  ignition  is  less  than 
that  of  pressed  lead  azide,  and  the  mean  time  required 
for  DPT  is  longer  than  that  of  the  lead  azide  but  shorter 
than  that  of  lead  styphnate. 
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(3)  The  critical  length  for  DDT  in  cast  2-MNT  is  approximately 
2 mm,  and  the  critical  time  required  for  the  full  develop- 
ment of  detonation  in  the  DDT  process  is  roughly  8 ixsec. 

C.  Modification  of  Macek's  Model  for  DDT  in  Cast  Explosives 

Early  work*°»*^  on  the  transition  from  deflagration  (slow  pressure- 
dependent  burning)  to  detonation  led  to  the  hypothesis  that  the  sequence 
of  events  in  the  transition  can  be  divided  into  three  parts:  a rapid 

pressure  increase  behind  the  burning  front  that  sends  compression  waves 
through  the  front  and  into  the  unreacted  explosive,  the  coalescence  of 
these  compression  waves  into  a shock  wave  in  the  unreacted  explosive 
some  distance  ahead  of  the  flame,  and  the  shock  init  ion  of  a detonation 
wave  in  the  unreacted  explosive.  This  sequence  explained  the  observations 
of  detonation  waves  in  gaseous  and  condensed  explosives  that  traveled 
back  toward  the  burning  front  from  the  point  of  initiation  of  detonation, 
Macek^^ > then  developed  a first-generation  experimental  and  theoretical 
approach  to  the  study  of  DDT  in  cast  secondary  explosives.  Because  this 
research  program  entails  the  study  of  DDT  in  cast  primary  explosives, 
an  analysis  of  Macek's  approach  and  the  extension  of  this  approach  to 
cast  primaries  were  undertaken.  This  section  describes  the  analysis  of 
Macek's  approach  and  the  derivation  of  an  improved  approach;  the  next 
section  describes  the  application  of  the  improved  approach  to  lead  azide. 

Macek's  approach  to  the  study  of  DDT  in  cast  secondaries  was  to 
make  experimental  measurements  of  the  pressure- time  history  in  the 
deflagration  region  and  of  the  time  and  position  where  detonation  first 
appeared.  The  apparatus  used  by  Macek  is  shown  in  Figure  9.  A strain 
gage  positioned  on  the  outside  of  the  tube  was  used  to  measure  the 
pressure-time  history  In  the  deflagrating  explosive  near  the  bridgewire. 
Ionization  probes  were  located  at  regular  Intervals  to  record  the  passage 
of  the  deflagration  wave  or  compression  waves  in  the  unreacted  explosive. 

The  measured  pressure- time  rise  was  of  the  form 
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where  P Is  the  initial  pressure,  t Is  the  time,  k Is  an  experimentally 
o 

determined  constant,  and  P Is  the  pressure  in  the  reaction  product  gases. 
This  pressure  was  assumed  to  be  equal  to  the  pressure  In  the  unreacted 
solid  explosive  which  obeyed  the  modified  Tait  compression  equation 


where  p Is  the  Initial  density,  c Is  the  Initial  sound  velocity, 
o o 

p Is  the  density  of  the  compressed  explosive,  and  n Is  a constant 
(generally  n|=  3).  The  propagation  of  the  flame  and  Its  associated 
high-pressure  gases  into  the  solid  explosive  creates  compression  waves. 

The  motion  of  these  compression  waves  through  the  unreacted  solid  was 
calculated  by  the  method  of  characteristics,  which  is  illustrated  in 
a distance-time  diagram  such  as  Figure  10.  According  to  the  Riemann 
analysis, the  quantity  u - (T  is  Invariant  along  a negative  characteristic 
(u  - c),  where  u is  the  particle  velocity  in  the  unreacted  solid, 


and 


Thus,  the  characteristics  in  Macek's  analysis  were  straight  lines.  As 
previously  pointed  outi  straight-line  characteristics  assume  that  no 
reaction  occurs  in  the  compression  region.  Because  the  initial  particle 
velocity  was  zero,  the  particle  velocity  was  given  by 

u(t)  = a(t)  - a = c(t)  - c (7) 

o o 

The  velocity  of  the  compression  waves  was  u(t)  + c(t)  and  the  position 
x(t)  of  the  gas-solid  boundary  was 

t 

x(t)  =1'  (c  - c )dt  (8) 

o 

When  the  velocities  of  the  compression  waves  were  calculated  at 
various  pressures,  and  drawn  in  Figure  10  as  characteristics  (u  + c), 
the  characteristics  intersected  at  a distance  of  about  12  cm  into  the 
unreacted  explosive.  Figure  11  shows  the  development  of  the  compression 
wave  in  pressure-distance  plots  at  various  times.  This  intersection  of 
the  characteristics  was  assumed  to  represent  the  formation  of  a shock 
wave  of  sufficient  magnitude  and  duration  to  cause  initiation  of  detonation 
in  that  region  of  the  sc' id  explosive.  Once  initiated,  the  detonation 
wave  traveled  through  the  unreacted  explosive  in  the  direction  of  flame 
propagation  (detonation)  and  also  toward  the  flame  front  (retonatlon) . 

Using  this  approach,  Macok  obtained  qualitative  agreement  between 
calculated  times  and  distances  required  for  DDT  and  observed  times 
(~  100  )isec)  and  distances  ('v.  12  cm)  for  two  cast  secondaries, 
diethylnitramlne  dinitrate  (DINA)  and  pentollte  50/50.*®  Macek  also 
attempted  to  explain  the  observed  exponential  pressure-time  relationship 
[equation  (1)1  behind  the  burning  front  on  the  basis  of  one-dimensional, 
adiabatic  flame  model.  This  flame  had  a burning-rate  law  of  the  form 
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FIGURE  11  DEVELOPMENT  OF  COMPRESSION  WAVE  FROM  DEFLAGRATION 
AS  TRANSPOSED  FROM  THE  CHARACTERISTICS  DIAGRAM 
IN  FIGURE  10 

The  left  boundary  of  the  shaded  area  represents  the  position  of  the 
defiagration  surface. 


where  R was  the  burning  rate,  0 was  a constant  equal  to  10  cm/sec-kbar, 
and  X was  a constant  equal  to  one.  This  type  of  burning- rate  law  has 
been  observed  fur  several  explosives  over  a range  of  ambient  pressures 
up  to  a kllobar  by  Andreev  and  Chulko.'*  Macek  normalized  the  pressure- 
time relationship  derived  for  this  burning-rate  law  to  equation  (1) 
through  an  adjustable  proportionality  constant  that  Included  the  surface 
area  of  the  burning  front. 

Although  Macek' s approach  gave  qualitative  agreement  with  DDT 
experiments  and  contained  many  of  the  physical  processes  Involved  In 
DDT,  It  has  several  weak  points  that  require  modification.  In  the 
experimental  setup,  Macek  used  strain  gages  placed  on  the  outside  of 
the  tube  to  measure  pressure  because  those  were  the  best  available  tools 
at  that  time.  Recent  advances  In  the  use  of  In-materlal  stress  gages 
make  these  crude  strain  gage  measurements  obsolete.  In  particular, 
multiple  Lagrange  stress  gages  have  been  successfully  used  In  reacting 
explosives*’  and  near  exploding  brtdgewlres  to  measure  detonator  output,*® 
The  recommended  use  of  these  stress  gages  to  measure  pressure- time 
histories  during  DDT  in  cast  2-MNT  Is  discussed  In  Section  IV,  Conclusions. 
In  addition  to  better  pressure  measurements,  the  deflagration  velocity 
measurements  must  be  much  more  accurate  than  In  Macek 's  model  to  determine 
the  actual  burning  rate  and  surface  area  of  burning  as  functions  of  time. 

On  the  theoretical  side  of  Macek* s approach,  three  problems  exist. 

By  assuming  that  the  pressures  in  the  burned  gases  and  in  the  unreacted 
explosives  were  equal,  Macek  neglected  the  pressure  drop  across  a 
deflagration  wave,  which  can  be  very  large  If  the  flame  approaches  a 
Chapman-Jouguet  (CJ)  deflagration.  Thus  the  pressure  in  the  unreacted 
explosive  ahead  of  the  flame  may  have  been  much  higher  than  the  pressures 
measured  near  the  brldgewlre  In  Macek *s  apparatus. 

The  second  problem  was  pointed  out  by  Jacobs,*®  who  noticed  that  the 
linear  pressure-burning  rate  relationship  derived  by  Macek  cannot  describe 


the  pressure-particle  velocity  buildup  at  pressures  above  one  kllobar. 

He  showed  that  Macek's  linear  pressure/burning  rate  relationship  can 
account  for  the  exponential  pressure/tlme  relationship  only  when  surface 
area  of  the  burning  front  Is  20,000  times  the  original  surface  area  of  the 
cast  explosive.  Even  the  Introduction  of  mechanical  means  of  Increasing 
the  surface  area  available  to  the  flame,  such  as  fracture  and  Induced 
porosity,  cannot  account  for  that  magnitude  of  surface  area  Increase. 
Therefore,  the  actual  burning  rate-pressure  relationship  In  accelerating 
deflagrations  that  cause  detonation  must  be  very  different  from  the  linear 
burning  rate^pressure  law  observed  for  explosives  burning  at  high  ambient 
pres^fares  and  low  ambient  temperatures.*® 


The  third  problem  with  Macek's  analysis  Is  that  the  momentum  Imparted 
to  the  burnt  gases  by  the  deflagration  wave  Is  neglected.  In  a closed 
tube,  this  momentum  participates  In  the  flame  acceleration  mechanism. 

Figure  12  shows  the  detonation  and  deflagration  branches  of  the  Hugonlot 
curves  In  the  pressure-volume  and  pressure-particle  velocity  planes.  In 
a detonation  wave  the  pressure  and  particle  velocity  Increase  as  the  volume 
decreases,  while  In  a deflagration  wave  the  pressure  and  particle  velocity 
decrease  as  the  volume  Increases.  When  the  deflagration  wave  Is  Ignited 
in  a closed  tube,  it  travels  up  the  tube  burning  the  explosive  and  creating 
a flow  of  burnt  gases  in  the  opposite  direction  (i.e.,  toward  the  real  wall). 
As  these  gases  contact  the  wall,  their  kinetic  energy  Is  converted  Into 
Internal  energy.  The  resulting  pressure  Increase  interacts  vlth  the 
deflagration  wave,  thereby  establishing  a mechanism  of  flame  acceleration. 
Because  the  accelerating  deflagration  wave  Is  creating  a comprci^sion  field 
in  the  unreacted  explosive,  a complicated  nonsteady  flow  d^elo^s.  This 
flow  consists  of  a compression  field  ahead  of  the  flame  In  w(' IcH.pi.rtlcles 


of  unreacted  explosive  have  a velocity  In  the  forward  direcilf  fnd 

accelerating  turbulent  flame  zone  In  which  the  unreacted  explorMe  burns, 

\ 

loses  its  forward  particle  velocity,  and  becomes  a burnt  gasj  moyture  with 


43 


> 

\ 

% 


¥ 


44 


a negative  particle  velocity;  and  a compression  zone  behind  the  flame  in 
which  the  burnt  gases  are  slowed  to  the  zero  particle  velocity  rear 
boundary  condition. 

This  complex  flow  is  illustrated  by  a series  of  particle  velocity- 
distance  diagrams  in  Figure  13.  This  type  of  flow  diagram  for  deflagration 
was  first  presented  by  Zeldovlch  and  Kompaneets, In  Figure  13,  the 
deflagration  wave  is  shown  as  a discontinuity  for  simplicity  but  it 
actually  has  a finite  width.  In  Figures  13a  and  13b,  the  positive 
particle  velocity  of  the  compressed,  unreacted  explosive  is  less  than  the 
negative  particle  velocity  imparted  by  deflagration,  so  a compression 
zone  is  required  to  bring  the  particle  velocity  to  zero  at  the  rear 
boundary.  Figure  13c  shows  the  unique  case  in  which  the  particle 
velocity  lost  in  deflagration  exactly  balances  the  particle  velocity 
previously  gained  in  compression.  At  later  times,  as  illustrated  in 
Figure  13d,  the  forward  momentum  is  not  totally  destroyed  by  deflagration, 
and  a rarefaction  wave  is  created  behind  the  flame  to  satisfy  the  rear 
boundary  condition.  However,  once  the  compression  zone  behind  the  flame 
no  longer  exists,  this  flame  acceleration  mechanism  also  vanishes.  The 
associated  compression  field  preceding  the  flame  must  then  be  strong 
enough  to  shock  initiate  the  unreacted  explosive  for  DDT  to  occur. 

Although  an  exact  solution  of  the  complex  nonstcady  flow  field 
resulting  from  deflagration  in  a closed  tube  is  extremely  difficult, 
two  analytical  models  of  flame  acceleration  that  provide  upper  and 
lower  deflagration  velocity  limits  have  been  formulated  by  considering 
the  nature  of  the  deflagration  process.  The  pressure-volume  diagram  in 
Figure  12  shows  that  both  the  detonation  and  deflagration  branches  of  the 
reaction  products'  Hugoniot  curve  have  CJ  states.  For  a detonation  wave, 
the  CJ  state  represents  the  only  stable  end  state  for  the  products,®* 
but  the  end  state  of  a deflagration  wave  may  be  represented  by  any 
point  on  the  weak  deflagration  branch  of  the  Hugoniot  curve,®*  which  is 
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FIGURE  13  PARTICi-E  VELOCITY-DISTANCE  DIAGRAMS  FOR  ACCELERATING 
DEFLAGRATION  WAVES 


bounded  by  points  A and  CJ'  in  Figure  13,  Point  A corresponds  to  a 
constant  pressure  deflagration  in  which  no  particle  velocity  is  produced 
and  point  CJ'  corresponds  to  the  maximum  possible  deflagration  and  particle 
velocities  and  pressure  drop.  Thus  the  final  state  of  a deflagration  is 
not  well  defined  as  it  is  In  a detonation,  but  steady  deflagration  waves 
in  tubes  tend  toward  the  CJ  state  at  long  times. 

Combining  the  flame  acceleration  mechanism  and  the  properties  of 
deflagration  waves  led  to  the  formulation  of  two  analytical  models  of 
flame  acceleration  that  bracket  the  actual  deflagration  process  in  a 
closed  tube.  One  model,  originally  used  by  Adams  and  Pack®®  for  DDT 
calculations  in  gases,  consists  of  the  compression  region  followed  by 
a deflagration  wave  whose  velocity  is  always  such  that  the  forward 
momentum  imparted  to  the  explosive  by  the  compression  waves  is  exactly 
canceled  out  by  the  loss  in  forward  momentum  caused  by  deflagration. 

Tha  Adams  and  Pack  model  does  not  require  a rarefaction  or  compression 
wave  to  satisfy  the  rear  boundary  condition,  and  all  of  its  states  are 
represented  by  Figure  13c.  Since  the  change  in  particle  velocity  with 
distance  behind  the  flame  is  zero,  the  change  in  pressure  at  the  rear 
boundary  with  time  must  also  be  zero.  Although  this  model  conserves 
momentum,  it  is  not  physically  correct  because  it  contains  no  flame 
acceleration  or  rear  wall  pressure- time  increase  mechanisms.  However, 
the  Adams  and  Pack  model  provides  lower  limits  of  the  d.ef3.agration 
velocity  and  pressure  drop  across  an  accelerating  deflagration  wave, 
because  any  slightly  faster  deflagration  produces  negative  particle 
velocity  and  thus  possesses  the  necessary  mechanisms. 

The  second  model  of  flame  acceleration  assumes  that  the  deflagration 
is  always  a CJ  deflagration  wi rh  the  largest  possible  velocity  and  pressure 
drop  associated  with  a certain  initial  pressure.  This  moael  definitely 
Includes  the  flame-acceleration  and  rear-wall  pressure-increase  mechanisms 
exhibiting  all  of  the  flow  situations  shown  in  Figure  13.  The  pressure- 
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particle  velocity  profiles  for  CJ  deflagrations  relating  to  Figure  13  are 
shown  ii.  Figure  14.  This  CJ  deflagration  model  was  discussed  for  DDT 
in  gases  hy  Troshin,^^  An  interesting  state  where  the  two  models  predict 
the  same  velocity  and  burnt  gas  pressure  is  represented  by  Figures  13c  and 
14c.  At  this  point  the  Adams  and  Pack  deflagration  becomes  a CJ  deflagration, 
and  the  CJ  deflagration  reaches  a velocity  for  which  the  momentum  loss 
. during  burning  exactly  cancels  the  forward  momentum  of  the  unreacted 
.explosive. 

For  the  CJ  deflagration,  this  state  also  corresponds  to  the  end  of 
the  flame  acceleration  mechanism  of  compression  waves  overtaking  the 
deflagration  wave.  Table  4 shows  the  calculated  properties  of  this  state 
behind  a shock  wave  in  cast  TNT  and  cast  pentolite  50/50.  The  resulting 
shock  Ttressures,  95  kbar  for  cast  TNT  and  129.5  kbar  for  cast  pentolite 
50/50,  are  too  high  to  indicate  a criterion  for  DDT,  because  pentolite, 
which  readily  undergoes  DDT,  has  a shock  sensitivity  of  about  20  kbar®® 

I 

iand  cast  TNT,  which  does  not  undergo  DDT,  can  be  initiated  by  a 60-kbar 

shock. Therefore,  some  other  criterion  for  flame  acceleration  determines 

'V 

whether  a cast  explosive  undergoes  DDT. 

These  two  analytical  models  of  flame  acceleration  were  used  to  make 
DD7  c ilct^fttions  for  cast  pentolite,  which  is  the  only  explosive  for  which 
a coraplei.T.'  the  required  data  was  available.  The  method-of- 

characteristics  pa»rt  i»f  Mpcek's  analysis  was  used  to  predict  the  time 
and  distance  to  shock  formation  ahead  of  the  flame.  The  pressure-time 
histories  used  were  those  of  Macek^^  up  to  one  kilobar 


P ^ (0.08  kbar)  e 


0.12t 


and  the  more  recent  data  of  Price  and  Wehner  above  one  kilobar 


P„  = <1  kbar)  e 
2 


0,08856(t-21.05) 
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Table  4 

CONDITIONS  REQUIRED  FOR  A CJ  DEFLAGRATION 
TO  A ZERO  PARTICLE  VELOCITY  END  STATE 
BEHIND  A SHOCK  WAVE 


1. 1 J .'V-  MWIMI 


i'j 


Cast  TNT 

Shock  Pressure  in  Explosive  (kbar)  95.0 

CJ  Deflagration  Velocity  (mrd/|Jsec)  3.116 

CJ  Pressure  after  Deflagration  (kbnr)  43.5 


Cast  Pentollte  50/50 
129.5 
3 , 504 
63.4 
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where  P Is  the  pressure  behind  the  deflagration  and  t Is  in  microseconds. 

A 

The  conditions  ahead  of  the  flame  were  determined  by  the  equations 
for  an  isentropic  compression  with  n = 3 In  the  modified  Talt  equation. 

The  equations  for  the  density  p^,  particle  velocity  u^,  characteristics 
(Uj^  + c^)  , and  Internal  energy  In  terms  of  the  pressure  p^  in  the 
compressed  unreacted  explosive  are 


Pq=  P^/l- 


(12) 


u = c 
1 o 


1 + 


1/3 


as) 


u,  + c,  -.=  c 
1 1 o 


1/3 


- 1 


(14) 


(15) 


The  conservation  equations  across  the  deflagration  wave  are 


(MASS)  p.  (W  - u ) = p^  (W  - u ) (16) 

X 12  2 

2 2 

(MOivIENTUM)  P + P (W  - u ) a P + p (W  - u ) (17) 

XX  ] 2 2 2 

Pi  (W  - u )^  P (W  - u )^ 

(ENERGY)  E + — + + Q = E + — + (18) 

1 P^  2 2 p^  2 

where  the  subscript  2 denotes  the  state  behind  the  deflagration  wave, 

W is  the  deflagration  velocity  and  Q is  the  heat  of  reaction.  The 
equation  of  state  for  reaction  products  is 
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and 


-K 

P p„  = constant 
2 2 


(19) 


, K 2 

E + P /p„  = •: — 

2 2 2 K - 1 p 

2 


(20) 


where  K is  the  adiabatic  expansion  coefficient  determined  at  the  CJ 
detonation  state  through  the  equation 


CJ 


K + 1 


= 2 (K  - 1)Q 


(21) 


where  D is  the  detonation  velocity. 

In  the  Adams  and  Pack  model  of  zero  particle  veloclt  behind  the 
deflagration,  u is  set  equal  to  zero  and  the  resulting  equation  for 

A 

the  deflagration  velocity  In  terms  of  state  1 Is 


(k-ihq.k, 

Ui  [ 1 p^  2 1 


(K-1) 


0 


(22) 


The  pressure  P and  density  p are  then  calculated  from  Equation  (16)  and 

A 4b 

(17)  with  u s=  0.  Table  5 lists  the  Input  data  on  pentollte  50/50 

A 

required  for  the  calculations  and  the  calculated  values  of  u^, 

W,  and  P for  various  values  of  P , the  pressure  in  the  compressed, 

A X 

unreacted  explosive.  As  shown  in  Table  5,  the  calculated  deflagration 

velocities  W are  roughly  300  times  faster  than  predicted  by  the  linear 

12 

pressure-burning  rate  relationship.  Equation  (9).  Gibson  and  Macek 

5 5 

actually  observed  deflagration  wave  velocities  of  10  - 2 x 10  cm/sec, 

so  these  calculated  velocities  appear  to  be  of  the  correct  magnitude. 

The  closeness  of  the  values  for  P and  P clearly  shows  that  the  flame 

12 

in  the  Adams  and  Pack  model  Is  a nearly  constant  pressure  deflagration 
(corresponding  to  point  A in  Figure  11}  for  pressures  below  20  kbar. 


Table  5 


INPUT  DATA  FOR  PENTOLITE  AND  RESULTS 
FOR  DEFLAGRATION  WAVES  IN  PENTOLITE  (ADAMS  AND  PACK  MODEL) 


DATA  ON  PENTOLITE 


Density  = 1.67  g/cm 

5 

C = 2.43  X 10  cm/sec 

° 5 

D = 7.47  X 10  cm/sec 


259  kbar 
2.598 

1.16  kcal/g 


(kbar) 

P (kbar) 
2 

(cm/ sec) 

W(cm/sec) 

+ c 

1_ 

(cm/sec) 

1 

0.9998 

3 

V4C  10 

3 

2.489  X 10 

2.479 

X 

5 

10 

2 

1.9984 

3 

4.031  X 10 

3 

5.023  X 10 

2.527 

X 

5 

10 

4 

3 988 

3 

9.481  X 10 

4 

1.0228  X 10 

2.620 

X 

5 

10 

7 

6.937 

4 

1.6151  X 10 

4 

1.8358  X 10 

2.753 

X 

5 

10 

10 

9.815 

4 

2,249  X 10 

4 

2.700  X 10 

2.880 

X 

5 

10 

20 

18.75 

4.171  X 10^ 

4 

5.700  X 10 

3.264 

X 

5 

10 

30 

26.36 

4 

5.864  X 10 

4 

8.860  X 10 

3.603 

X 

5 

10 

40 

32.56 

4 

7.385  X 10 

5 

1,201  X 10 

3.907 

X 

5 

10 

70 

44.47 

5 

1.1244  X 10 

5 

2.054  X 10 

4.678 

X 

5 

10 

100 

50.28 

5 

1.4410  X 10 

5 

2.738  X 10 

5.311 

X 

5 

10 

123.8 

52.62 

5 

1.6597  X 10 

5 

3,186  X 10 

5.749 

X 

5 

10 

(CJ  state) 
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The  only  other  quantity  required  for  a distance-time  plot  for  DDT 


similar  to  that  In  Figure  8 for  pentollte  using  this  model  Is  the 
position  of  the  flame  front  as  a function  of  time,  x(t),  which  Is  derived 
from 


x(t)  = J W(t)  dt 


(23) 


where  W(t)  Is  determined  from  the  dependence  of  W on  P In  Table  4 and 

2 

the  .dependence  of  P on  time  from  Equations  (10)  and  (11).  Figure  13 

2 

shows  the  resulting  distance-time  plot  for  pentollte  using  the  Adams 
and  Pack  model  of  deflagration  behind  an  Isentroplc  compression.  Comparing 
Figure  13  with  Figure  8,  the  Adams  and  Pack  model  predicts  a much  more 
rapid  deflagration  and  thus  a greater  penetration  of  the  flame  Into  the 
explosive.  However,  the  time  and  distance  required  for  shock  formation 
aie  approximately  the  same  In  both  diagrams. 


These  DDT  predictions  must  also  be  compared  with  those  for  a CJ 
deflagration  following  an  isentroplc  compression.  At  the  CJ  state, 


= ^2  + c^,  or 


(W  - u )^  = K ~ 
" ^2 


(24) 


Substituting  Equation  (24)  into  the  conservation  Equations  (16)  - (18) 
and  solving  for  P^  yields  the  following  relationship 


P..  = (K-1)  p, 
£,  ' 


I (2Q  + 2E  + P /p  )P 

, (Q  + E ) 1 i = — = =• 

1 1 ■%/  2 2 
f P (K~1)(Q+E> 


and  solving  for  p gives 
2 


-1 


(K-1)(Q+E 


5+E^  + P^/P^)  r f K^(2Q  + 2E  + P /p  )P 

"kP  ( ^ ^ W ^ ” 2 ~2 

1 L ▼ -1)(Q+E^  + P^/P^) 


(25) 


(26'> 
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FIGURE  15  DISTANCE-TIME  DIAGRAM  FOR  SHOCK  FORMATION  IN  CAST  PENTOLITE  50/50  USING  THE  ADAMS 
AND  PACK  MODEL  OF  DEFLAGRATION 


The  particle  velocity  behind  the  deflagration  and  the  deflagration 
velocity  W are  then  calculated  from 


W - 


and 


(27) 


(28) 


Table  8 lists  the  calculated  states  for  a CJ  deflagration  in  Isentroplcally 
compressed  pentollte  at  various  pressures.  Figure  16  Is  the  distance- 
time plot  for  DDT  for  the  CJ  deflagration  model.  The  pressures  behind 
the  CJ  deflagration  wave  are  much  lower  than  the  corresponding  values  for 
the  Adams  and  Pack  model  at  low  initial  pressures  P^,  and  the  CJ 
deflagration  velocities  are  nearly  twice  as  large  as  those  In  Table  5. 

In  terms  of  the  time  and  distance  to  DDT,  the  CJ  deflagration  front  moves 
a greater  distance  Into  the  explosive,  and  the  shock  formation  occurs  In 
a shorter  time  (about  80  fisec  In  Figure  16  compared  with  90  jUsec  In 
Figure  15),  whereas  the  distance  to  shock  formation  Increases  (about 
14  cm  In  Figure  16  compared  with  about  13  cm  In  Figure  15). 

Since  these  two  analytical  models  represent  upper  and  lower  limits 
to  the  deflagration  rate,  the  critical  time  and  distance  for  DDT  In  cast 
pentollte  50/50  should  be  bracketed  by  these  results  If  the  Isentroplc 
compression  of  cast  pentollte  50/50  obeys  equation  (2)  with  n = 3,  More 
fully  Instrumented  DDT  experiments  on  pentollte  are  required  to  test  the 
validity  of  these  results  but  the  distance,  times,  pressures,  and 
velocities  to  be  measured  are  defined  by  the  following  results: 

(1)  Rear  wall  pressure  history  measurements  of  20-40 
kbar  in  the  first  50  fisec  after  ignition. 

(2)  Deflagration  velocities  of  10^-10®  cm/sec  In  the 
first  2 cm  of  explosive. 
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Table  6 


CJ  DEFLAGRATION  WAVE  CALCULATIONS 
FOR  ISENTROPICALLY  COMPRESSED  PENTOLITE 


Pj^(kbar) 

Pg  (kbar) 

u (cm/sec) 
2 

W(cm/sec) 

1 

0.2799 

-2.044  X 10® 

3 

4.500  X 10 

2 

0.5637 

-2.019  X 10® 

8.910  X 10® 

4 

1.1426 

-1.958  X 10® 

1.7503  X 10 

7 

2.038 

-1.870  X 10® 

2.986  X 10^ 

10 

2.966 

-1.782  X 10® 

4.170  X 10^ 

20 

6.261 

-1.5353  X 10® 

7.767  X 10*^ 

30 

9.835 

-1.3208  X 10® 

1.0962  X 10 

40 

13.65 

-1.1304  X 10® 

1.3862  X 10 

70 

26.32 

-6.539  X 10“* 

2.130  X 10® 

100 

38.54 

-2.247  X 10^ 

2.719  X 10® 

123.8 

52.62 

0 

3.186  X 10® 
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DEFLAGRATION  WAVE 


(3)  Shock  formation  measurements  at  13-14  cm  into  the 
explosive  at  times  of  80-90  ijsec  after  Ignition, 

DDT  calculations  using  the  Adams  and  Pack  and  CJ  deflagration  models 
In  the  Macek  approach  yield  much  more  realistic  results  than  the  original 
burning  rate  model.  In  the  next  section,  these  models  are  used  to  make 
DDT  calculations  on  lead  azide,  the  only  primary  for  which  some  meager 
data  exist, 

D.  DDT  and  STD  Calculations  for  Lead  Azide 

The  processes  of  SDT^®  and  DDT*®  have  been  quantitatively  studied  to 
a limited  extent  In  lead  azide.  Chaudhrl  and  Field*®  observed  very  fast 
deflagration  waves  In  hot-wire  Initiation  tests  on  lead  azide  crystals 
of  thicknesses  below  the  critical  thickness  for  DDT.  By  photographing 
the  deflagrating  crystals  under  water,  they  measured  the  velocities  and 
pressures  of  the  shock  waves  produced  In  the  surrounding  water  by  the 
rapid  deflagration  of  lead  azide.  These  pressure  data,  combined  with 
Chaudhrl  and  Field's  framing  camera  records  of  the  Ignition  process, 
provide  a rough  estimate  of  the  pressure-time  history  In  the  deflagration 
region.  Chaudhrl 's  measurements®®  of  the  sound  velocity  and  shock  Initia- 
tion of  lead  azide  crystals  and  Davies'  measurements*®  of  the  unreacted 
Hugonlot  and  sound  velocity  In  lead  azide  pressed  to  a density  of 
3.4  g/cm®  provided  data  on  the  properties  of  lead  azide  In  the  compression 
zone  ahead  of  the  deflagration  wave.  Usli^g  these  data  In  the  two  models 
of  deflagration  discussed  In  the  previous  section,  we  calculated  the 
times  and  distances  for  DDT  In  lead  azide  crystals  and  pressed  lead 
azide  and  for  SDT  In  pressed  azide  and  compared  them  with  the  available 
experimental  results  as  a preliminary  test  of  the  DDT  models. 

To  calculate  the  time  and  distance  required  for  DDT  In  lead  azide 
crystals,  we  first  related  the  pressure  data  obtained  by  Chaudhrl  and 
Field*®  to  the  time  after  Ignition.  The  actual  pressures  measured  were 

59 


those  of  the  shock  waves  produced  in  water  as  the  lead  azide  deflagration 
front  reached  the  end  of  the  crystal.  Using  the  acoustic  approximation 
developed  by  Deal®’  for  a detonation  wave  propagating  into  an  inert 
material,  th'^^  pressure  in  the  lead  azide,  is  given  by 


P 

W 


+ 


P U \ 
LA  LA  \ 

P u I 
W / 


(29) 


where  P is  the  measured  water  pressure,  p is  the  density  of  the  water, 

w w 

U is  the  measured  shock  velocity  in  the  water,  p is  the  density  of 
w 3 LA 

or-lead  azide  crystals  (4.60  g/cm  ),  and  is  the  measured  velocity  of 

propagation  of  thexeactlon  front  in  the  lead  azide.  Tabiu  7 lists  the 

measured  values  and  the  calculated  pressures  in  the  lead  azide  crystals. 

Since  P is  the  pressure  at  the  front  of  the  wave,  it  corresponds  to 
LA 

P^,  the  pressure  in  the  unreacte.^  explosive.  These  pressures  were  related 
to  the  time  after  ignition  by  examination  of  Chaudhri  and  Field's  framing 
camera  records  of  these  reactions.  I.i  these  records  the  framing  rate  was 
0.7  M<sec  between  frames  and  the  transition  from  ignition  to  full  develop- 
ment of  the  steady  deflagration  velocity  for  a certain  crystal  thickness 

took  less  than  one  frame.  Thus  the  calculated  pressures  P . were  assumed 

LA 

to  be  attained  in  less  then  0.7  M>sec.  The  largest  pressure  P (12.85  kbar) 
was  assumed  to  be  dt  3loped  in  exactly  0.7  ^isec  and  the  c r two  pressures 
were  assumed  to  occur  at  times  of  0.7  p-sec  minus  the  difference  in  crystal 
thicknesses  divided  by  the  average  deflagration  velocity  (2.7  mm/pisec)  in 
these  experiments,  The  resulting  pressure-time  curve  was  of  the  fora 


. 6.94t 

P . = P,  = (0.1  kbar)e 
LA  1 


(30) 


where  t is  in  microseconds. 

The  particle  velocity  u^  in  the  unreacted  lead  azide  was  related 
to  P^  using  the  relationship 


- : . . r , U-  , 

iUfttiiiiiajliilliiiliKiiii 
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p,  = p C u 

1 o o 1 


(31) 


which  was  found  by  Davies  et  al.‘®  to  hold  for  pressed  lead  azide  at 

pressures  below  10  kbar.  Chaudhri^®  reported  a propagation  rate  for  a 

low  amplitude  longitudinal  stress  wave  in  lead  azide  crystals  of 

2.15  ± 0.3  mm//isec.  This  value  was  taken  for  c In  these  calculations. 

o 

Actually,  if  the  sound  velocity  is  2.15  mm//isec,  the  waves  that  Chaudhri 
and  Field refer  to  as  fast  deflagration  waves  are  supersonic  and  therefore 
are  some  sort  of  low-velocity  detonation  waves  (LVD)  that  cannot  develop 
into  high-velocity  detonation  (HVD)  because  of  rapid  pressure  decay  at 
the  edges  of  the  crystal.  A phenomenon  like  LVD  has  been  observed  for 
lead  azide  pressed  to  high  densities  (3.75  g/cm®)  by  Leopold.®® 


The  remaining  features  of  the  DDT  calculations  are  identical  to  those 

described  in  the  last  section.  Table  8 lists  the  input  data  for  crystalline 

lead  azide  and  the  results  of  calculations  at  various  initial  pressures 

for  the  two  deflagration  models,  the  Adams  and  Pack  model  and  the  CJ 

deflagration  model.  The  calculated  deflagration  velocities  for  both 

models  are  an  order  of  magnitude  lower  than  those  observed  by  Chaudhri 

and  Field,  while  the  velocities  of  C characteristics  (u  + c ) are  of 

5 11 
that  magnitude  (~2.5  x 10  cm/sec).  Figures  17  and  18  show  the  distance- 

time  diagrams  for  DDT  in  o'-lead  crystals  predicted  by  the  Adams  and  Pack 

and  the  CJ  deflagration  models,  respectively.  The  calculated  distances 

for  DDT  are  approximately  2 mm.  This  agrees  with  the  observed  behavior 

of  lead  azide  crystals  of  varying  thickness.  Crystals  less  than  2 mm 

thick  cannot  undergo  the  DDT  that  thicker  crystals  readily  make.  The 

calculated  times  to  DDT  are  approximately  2 ;j,sec  and  thus  agree  with 

the  preliminary  result.,  of  Slagg,®®  who  has  measured  times  for  DDT  of 

3±2  M^sec  for  crystals  of  lead  azide. 

The  DDT  calculations  for  crystalline  lead  azide  exhibit  order-of- 
magnitude  agreement  with  the  experimental  results  and  this  is  all  that 
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Table  7 

PRESSURES  AND  WAVE  VELOCITIES  FROM  CHAUDHRI  AND  FIELD'S 
EXPERIMENTS  ON  a-LEAD  AZIDE  CRYSTALS 


Crystal 

Thickness 

(am) 

Lead  Azide 
Deflagration 
Velocity 
(nun/asec) 

Water  Shock 
Velocity 
(mm/asec) 

Water  Shock 
Pressure 
(kbar) 

Calculated  Lead 
Azide  Pressure 
(kbar) 

310 

2.50 

1.67 

1.3 

5.14 

560 

2.90 

1.72 

2.0 

8.76 

670 

2.62 

1.90 

3.5 

12.85 
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can  be  expected  at  the  present  time.  The  exact  nature  of  the  rapid 
deflagration  or  LVD  waves  in  lead  azide  crystals  is  unknown.  The 
framing  camera  records  of  the  process  indicate  that  some  reaction 
occurs  at  the  front  of  the  wave.  The  general  approach  to  DDT  appears 
to  be  correct  but  the  effect  of  reaction  on  the  characteristics  used  in 
the  shock  formation  part  of  the  model  needs  to  be  considered. 

The  recent  SDT  experiments  on  pressed  lead  azide  by  Davies  et  al,^® 
pioneered  the  use  of  stress  gages  in  thin  flyer  plate  impacts  of  primary, 
explosives  and  also  yielded  data  on  the  unreacted  Hugoniot  and  on  the 
time  and  distance  for  SDT.  For  lead  azide  pressed  to  a density  of  3.4  g/cm®, 
Davies  found  that  stress  waves  propagate  at  1.23  mm/;isec,  and  the  equation 
for  the  Hugoniot  curve  in  this  case  is  obtained  as 

p^  =41.8  (32) 

by  substituting  these  values  of  c and  o into  equation  31.  For  both 

o o 

long  pulses  (3,5  ^sec)  and  short  pulses  (0,1  jUsec),  in  the  6-9  kbar 
range,  Davies  observed  that  the  transition  to  detonation  took  about 
2 Msec  and  occurred  about  1-2  ram  into  the  sample  of  pressed  lead  azide. 

To  determine  whether  the  theoretical  approach  to  DDT  developed  in 
this  report  would  predict  similar  times  and  distances,  we  used  this 
unreacted  Hugoniot  for  pressed  lead  azide  and  the  pressure-time  relation- 
ship of  equation  (30)  in  the  CJ  deflagration  model  to  calculate  the 
distance- time  diagram  shown  in  Figure  19,  This  calculation  predicts 
a shorter  time  (1.2  /isec)  and  distance  (0,9  mm)  for  pressed  lead  azide 
than  for  crystalline  lead  azide.  This  is  mainly  due  to  the  lower  sound 
velocity  in  pressed  lead  azide,  which  allows  the  compression  waves  to 
overtake  each  other  more  quickly.  To  compare  a prediction  of  the  CJ 
deflagration  model  to  Davies'  results,  the  SDT  calculation  shown  in 
Figure  20  was  made  for  an  8-ktaar  shock  that  causes  the  deflagration  to 
begin  about  0.6  usee  later  in  agreement  with  equation  (30).  This 
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INPUT  DATA  FOR  a-LEAD  AZIDE  CRYSTALS  AND  DEFLAGRATION  VELOCITY  CALCULATIONS 


EAD  AZIDE  CRYSTALS  USING  ADAMS  AND  PACK 


FIGURE  18  DISTANCE-TIME  DIAGRAM  FOR  DOT  IN  a-LEAD  AZIDE  CRYSTALS  USING  A CJ  DEFLAGRATION  MODEL 


DEFLAGRATION  MODEL 


FIGURE  20  DISTANCE-TIME  DIAGRAM  FOR  SDT  IN  TRESSED  LEAD  AZIDE 


calculation  resulted  in  the  formation  of  a second  shock  wave  due  to  the 
crossing  of  the  characteristics  ahead  of  the  flame.  This  second  shock 
wave  may  cause  detonation  at  1.3  fisec  and  1.1  mm,  or  detonation  may 
result  when  the  second  shock  overtakes  the  initial  shock  at  2.3  i|Sec  and 
2.8  mm.  Some  of  Davies'*®  gage  records  did  appear  to  detect  two  finite 
strength  shocks  at  1 or  2 mm  into  the  explosive  in  tests  where  detonation 
occurred  at  3 or  4 mm.  Again  the  order  of  magnitude  agreement  between 
the  SDT  approach  and  Davies'*®  results  is  as  good  as  can  be  expected  at  the 
present  time. 

From  these  DDT  and  SDT  calculations  on  pentolite  and  lead  azide,  it 
is  apparent  that  the  experimental  data  on  pressure-time  histories  in  the 
deflagration  region,  on  deflagration  velocity,  and  on  shock  formation 
ahead  of  the  flame  are  scarce  and  poor.  Thus  a good  test  of  the  modified 
Macek  approach  to  DDT  cannot  be  made  on  primaries  or  secondaries.  Well- 
instrumented  DDT  experiments  on  primary  explosives  are  critically  needed 
to  provide  a first  step  toward  an  understanding  of  DDT  in  primaries. 
Experiments  using  stress  gages  and  wave  velocity  measurement  techniques 
during  DDT  in  cast  2-MNT  and  other  primaries  are  suitable  for  providing 
precise  pressure  and  velocity  data  to  use  in  further  development  of  the 
theory. 


IV  CONCLUSIONS 


During  the  last  period  of  the  contract  series,  we  determined  the 
threshold  voltages  for  initiation  of  cast  1-MNT  and  2-MNT.  The  time  for 
DDT  as  a function  of  charge  length  was  measured  for  cast  2-MNT,  Two 
analytical  models  of  deflagration  were  Incorporated  Into  Macek's  theoreti- 
cal description  of  DDT,  resulting  in  Improved  calculations  for  DDT  In 
cast  pentollte.  These  models  were  also  tested  to  a limited  extent  on 
the  meager  existing  data  on  lead  azide.  In  analyzing  the  experimental 
work  on  DDT  In  cast  secondaries,  the  lack  of  precise  pressure  and  velocity 
measurements  was  apparent. 

In  future  research,  precise  measurements  of  pressure  and  wave 
velocities  should  be  made  In  the  critical  regions  of  cast  primaries  where 
DDT  develops.  These  measurements  will  give  the  first  detailed  results 
of  this  type  for  cast  primaries  or  secondaries.  These  data  will  furnish 
the  necessary  input  for  a realistic  test  of  the  theoretical  DDT  approach. 

The  precise  pressure  measurements  will  entail  the  use  of  Lagrange 
stress  gages  in  an  apparatus  similar  to  the  tlme-for-DDT  measurement 
apparatus.  Three  recent  applications  of  these  gages  in  explosives  re- 
search imply  that  their  use  in  DDT  of  primaries  will  be  successful.  As 
previously  discussed,  Davies  et  al.^®  successfully  used  manganln  and 
carbon  stress  gages  in  flyer  plate  impact  SDT  studies  of  pressed  lead 
azide.  In  recent  gage  development  work  at  SRI,  manganln  stress  gages 
were  used  to  measure  detonator  output  at  distances  of  a few  millimeters 
from  an  exploding  bridgewire. In  DDT  experiments  only  a heated  bridge- 
wire  Is  used,  so  the  brldgewlre  noise  level  should  be  even  less  of  a 
problem.  Mangani.  stress  gages  have  also  survived  reactive  shock  pres- 
sures of  60-80  kbar  in  shock  initiation  studies  of  cast  TNT.®'^  These 
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pressures:  are  higher  than  the  range  of  pressures  to  be  measured  during 
the  DDT  process,  because  the  shock  Initiation  threshold  of  a primary 
should  be  about  10-20  kllobars.  Single  crystals  of  lead  azide  have  been 
Initiated  by  thin  flyer  plates  accelerated  to  velocities  corresponding 
to  10-kbar  pressures  In  the  crystal.*® 

The  stress  gage  experiments  will  measure  stress  histories  In  the 
deflagration  zone  near  the  bridgewlre,  In  the  transition  zone  where 
shock  formation  occurs,  and  In  the  detonation  region.  Although  pressure 
measurements  In  the  region  of  the  CJ  state  are  not  feasible  at  this  time, 
pressures  of  60-80  kllobars  WA„ld  certainly  Indicate  that  detonation 
will  occur.  The  deflagration  velocities  near  the  bridgewlre  will  also 
be  monitored  to  determine  whether  the  high  velocities  predicted  by 
our  analytical  models  are  reached  and,  if  these  velocities  are  reached, 
which  of  the  two  models  (Adams-Pack  or  CJ  deflagration)  predicts  the 
velocity  more  accurately.  The  experiments  will  be  carried  out  on  cast 
2-MNT  and  cast  1-MNT  to  determine  the  difference  in  pressure  histories 
that  makes  2-MNT  a much  better  primary  than  1-MNT,  They  may  also  be 
carried  out  on  other  available  primaries,  such  as  pressed  lead  azide. 

The  pressure  and  deflagration  velocity  histories  will  also  provide  the 
input  data  for  a valid  test  of  the  theoretical  description  of  DDT  in 
cast  explosives. 

Once  the  DDT  aspect  of  primary  explosives  Is  quantified,  the  main 
emphasis  of  the  research  program  can  be  placed  on  determining  the  structural 
properties  that  control  the  pressure-time  and  deflagration  velocity-time 
histories.  An  understanding  of  these  histories  from  a chemical  kinetics 
point  of  view  is  the  next  stop  In  understanding  the  action  of  primary 
explosives.  This  step  is  a difficult  one  to  take,  however.  One  possible 
approach,  for  primaries  like  lead  azide  that  decompose  Into  simple  reaction 
product  mixtures  (Pb  and  N^),  would  be  to  calculate  the  rate  of  unl- 
molecular  decomposition  from  the  measured  pressure  history  behind  the 
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deflagration  wave.  The  measured  deflagration  velocity  history  and  pressure 
history  In  the  unreacted  explosive  would  determine  the  pressure  drop  across 
the  deflagration  zone.  A realistic  equation  of  state  for  the  compressed 
unreacted  explosive  ahead  of  the  flame  front  would  complete  the  thermo- 
dynamic description  of  the  primary  explosive  as  a function  of  time  and 
distance.  This  thermodynamic  description  would  be  Incorporated  Into  a 
chemical  kinetic  model  of  unlmolecular  decomposition  that  Is  consistent 
with  the  experimental  results.  A similar  derivation  of  a chemical  kinetic 
model  for  decomposition  of  a more  complex  molecule  containing  C,  H,  0, 
and  N atoms,  such  as  2-MNT,  would  require  assumptions  about  the  sequence 
of  Intermediate  reactions  Involved  in  the  decomposition  process.  It  would 
also  be  Important  to  explain  the  measured  deflagration  velocity  history 
using  the  theory  of  conductive  flames  and  to  calculate  the  energy  release 
rate  within  the  deflagration  zone. 

The  final  step  In  developing  a useful  correlation  between  chemical 
structure  and  primary  explosive  behavior  Is  to  relt te  chemical  kinetics 
to  chemical  structure.  This  type  of  relationship  between  kinetics  and 
structure  has  been  developed  on  a qualitative  basis  for  many  chemical 
systems,  Including  liquid  secondary  explosives.^®  These  qualitative 
relationships  have  proven  useful  for  certain  applications  of  secondary 
explosives.  A logical  approach  for  primary  explosives  is  to  develop 
qualitative  relationships  between  the  chemical  kinetic  models  that  are 
consistent  with  DDT  experimental  results  and  the  chemical  structures  of 
various  types  of  primaries.  These  qualitative  relationships  may  provide 
sufficient  guidelines  for  producing  better  primaries.  However,  because 
of  the  complex  and  varied  nature  of  primary  explosives,  it  may  be  necessary 
to  generate  a great  deal  of  experimental  DDT  and  kinetic  data.  These  data 
would  serve  as  the  required  base  for  formulating  detailed  quantitative 
relationships  between  chemical  kinetics  and  structure,  which  include  con- 
tributions from  all  groups  present  In  the  molecule,  that  may  be  necessary 
to  explain  the  effects  of  structure  on  primary  explosive  performance. 
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Appendix  A 

SYNTHESES  OF  1-  AND  2-METHYL- 5 -NITROTETRAZOLES 

The  procedures  described  below  were  used  to  synthesize  1-  and  2- 
methyl-5-nitrotetrazole  for  use  in  the  measurements  of  deflagration  to 
detonation  transitions  of  these  cast  primaries. 

l-Methyl-5-Nitrotetrazole 

A mixture  of  50  g of  aminotetrazole  monohydrate,  19.1  g of  sodium 
hydroxide,  and  97  ml  of  water  was  heated  to  92-95°C,  and  31,7  g of 
dimethyl  sulfate  was  added  dropwise  over  40  minutes.  The  temperature  of 
the  reaction  mixture  was  held  at  95°C  for  an  additional  hour  and  then 
cooled  to  5°C  overnight.  The  product,  l-methyl-5-aminotetrazole,  was 
removed  by  filtration,  washed  with  cold  water,  and  dried  to  yield  24.2  g. 
This  synthesis  was  developed  by  Henry  and  Finnegan.’ 

A mixture  of  5 g of  l-methyl-5-aminotetrazole  from  above,  17  ml  of 
96%  sulfuric  acid,  and  300  ml  of  water  was  added  dropwise  over  1 hour 
to  a solution  of  53  g of  sodium  nitrite  in  500  ml  of  water  at  45°C. 

The  reaction  mixture  was  stirred  an  additional  hour  at  45°C  and  then 
cooled  to  25°C.  The  reaction  was  extracted  with  three  50-ml  portions  of 
methylene  chloride.  The  methylene  chloride  extract  was  washed  once  with 
5%  sodium  bicarbonate  and  once  with  water.  It  was  dried  over  magnesium 
sulfate  and  evaporated,  leaving  1.5  g of  yellow  oil.  This  was  combined 
with  6.0  g of  oil  from  another  reaction  and  dissolved  in  a minimum 
volume  of  diethyl  ether.  The  ether  solution  was  passed  through  a neutral 
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R.  Henry  and  W.  Finnegan,  J.  Amer.  Chem.  Soc.  76,  923  (1954). 


d V-:-. 


77 


m 


ip,^JUJl  liMWi^pi^ff  Wl 


alumina  column.  The  product  from  the  column  was  recrystallized  twice 

from  ether  to  yield  3.5  g of  l-methyl-5-nitrotetrazole,  mp  52°C.  No 

trace  of  nitrosamlne  could  be  detected  by  infrared  analysis.  The  product,  1- 

1- methyl-5-nitrotetrazole,  was  stored  in  methylene  chloride  solution  at 
0°C  in  the  dark.  This  synthesis  procedure  was  developed  by  Bagal  et  al.® 

2- Methyl-5-Nltrotetrazole 

A mixture  of  50  g of  aminotetrazole  monohydrate,  15  ml  of  96%  HjSO^, 
and  970  ml  of  water  was  added  over  45  min  to  a mixture  of  77.7  g of 
sodium  nitrite,  48.5  g of  copper  sulfate,  and  730  ml  of  water.  The 
temperature  was  allowed  to  rise  slowly  from  0°  to  25°C  during  the  addition 
period.  The  reaction  mixture  was  stirred  an  additional  2 hours  until 
nitrogen  evolution  had  stopped.  Concentrated  sulfuric  acid,  44  ml, 
was  added  and  the  mixture  was  cooled  to  5°C.  The  product,  copper  5- 
nitrotetrazole,  was  removed  by  filtration  and  washed  with  water.  This 
procedure  was  developed  by  von  Herz®®'  and  has  been  modified  by  Gilligan 
and  Kamlet.^'^ 

The  copper  salt  from  above  was  dissolved  in  boiling  water,  and  a 
sodium  hydroxide  solution  was  added  slowly  until  the  pH  was  ~ 9.  After 
boiling  for  30  min  more,  the  mixture  was  cooled  to  25°C  and  filtered  to 
remove  the  copper  oxide.  The  filtrate  was  concentrated  and  cooled  to 
5°C.  The  crystallized  sodium  salt  was  removed  by  filtration  and  air 
dried  to  yield  35  g.  This  procedure  was  developed  by  von  Herz.®^ 


*L.  I,  Bagal,  et  al.,  Khim,  Geterotslkl. , Soedin.  (1970),  p.  259; 

Chem.  Abstracts  T^,  111383h. 

von  Herz,  U.S,  Patent  2,066,954  (1937);  W.'.  H.  Gilligan  and  M.  J. 
Kamlet,  "Synthesis  of  Mercuric  5-Nitrotetrazole,"  NSWC/WOL/TR  76-146 
(6  December  1976). 
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A mixture  of  18  g of  sodium  5-nitrotetrazole  tetrahydrate  from 
above,  400  ml  of  water,  160  ml  of  acetone,  and  20  g of  methyl  iodide  was 
refluxed  for  3.5  hr.  More  methyl  iodide,  6 g,  was  added  after  the  first 
2 r.  The  acetone  was  removed  by  distillation,  and  400  m3  of  benzene 
was  added.  The  mixture  was  then  washed  with  100  ml  of  2%  sodium  hydroxide 
in  brine.  After  two  additional  washings  with  25-ml  portions  of  brine,  the 
solution  was  dried  over  magnesium  sulfate  and  evaporated  to  dryness.  The 
residue  was  dissolved  in  a mixture  of  200  ml  of  benzene  and  60  ml  of 
petroleum  ether  and  cooled  to  0°C  for  24  hr.  The  product  was  removed 
by  filtration  and  dried  under  vacuum.  Yield  was  4.8  g,  mp  85°C.  The 
product,  2-methyl-5-nitrotetrazole,  was  stored  in  methylene  chloride 
solution  at  0°C  in  the  dark.  This  procedure  was  developed  by  Henry  and 
Finnegan. ^ 
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R.  Henry  and  W.  Finnegan,  J.  Amer.  Chem.  Soc.  76,  923  (1954). 
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Appendix  B 

ESTIMATION  OF  HEATS  OF  FORMATION  AND  DENSITIES  OF 
1-  and  2-METHYL-5-NITROTETRAZOLE 

The  use  of  additivity  methods  to  estimate  the  heats  of  formation  at 
25°C  of  l-methyl-5-nitrotetrazole  (1-MNT),  2-methyl-5-nitrotetrazole 
(2-MNT),  and  for  comparison,  5-azidotetrazole  (5AZT)  is  described  below. 
All  heats  are  in  kcal/mole. 

The  heats  of  formation  of  1-MNT  and  2-MNT  were  obtained  by  estimating 
the  effect  of  adding  methyl  groups  to  5-nitrotetrazole  (TNO^)  at  the  1- 
and  2-positlons.  The  heat  of  formation  of  TNOg  was  obtained  by  assuming 
that  the  difference  in  the  heats  of  formation  of  TNOj  and  nitrobenzene 
(PhNOj)  is  the  same  as  the  average  difference  in  heats  of  formation  of 
tetrazole  compounds  (TX)  substituted  at  the  5-position  and  phenyl 
compounds  wi+-'  che  same  substituents  (PhX),  That  is,  the  basic  assumption 
is  that 

AH^(TNOa)  - ^iP^(PhNOj)  = - [fljVhX)  (B-1) 

f f f f 

The  values  for  tetrazole  compounds  were  obtained  from  the  work  of 

Domalski;^  the  values  for  phenyl  compounds  are  from  Cox  and  Pilcher.® 

The  average  value  of  rZlH^(TX)  - ffi  (PhX)  1 was  calculated  as  shown  in 

f F ' 

Table  D-1. 


^ E.  S.  Domalskl,  J.  Phys.  Chem.  Ref.  Data  I,  221  (1972). 

® J.  D.  Cox  and  G.  Pilcher,  Thermochemistry  of  Organic  and  Organometallic 
Compounds  (Academic  Press,  New  York,  1970). 
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Table  B-1 

COMPARISON  OF  THE  HEATS  OF  FORMATION  OF 
DERIVATIVES  OF  TETRAZOLE  AND  BENZENE 

Tetrazole  Derivatives  Benzene  Derivatives 
TX  PhX 


X _ 

(all  solid) 

(all  ideal  gas) 

A(TX  - PhX) 

H 

56,7 

19.8 

36.9 

(M 

1.5 

-23.0 

24.5 

NHs 

49,7 

20.8 

28.9 

CN 

96.1 

51.5 

44.6 

CH3O 

16.6 

-17,3 

33.9 

Average  34  ± 10 


Inserting  the  average  value  for  AH^(TX,  solid)  minus  ^^(PhX,  ideal  gas) 
into  equation  (B-1)  gave 

^^(TNOj,  solid)  - ^^(PhNOa,  ideal  gas)  = 34  (B-2) 

A value  of  16,9  kcal/mol  for  AH^CPhNOj),  from  previous  work,®  was  used 
to  give; 

^^(TNOj,  solid)  = 34  + 16,9  = 51  (B-3) 

The  effect  of  adding  methyl  groups  to  TNOj  at  the  1-  and  2-positions 
to  form  IMNT  and  2MNT  was  estimated  by  comparing  5-aminotetrazole  (AT)  with 
1-  and  2-methyl-5-aminotetrazole  (LMAT  and  2MAT),  That  is,  we  assumed 
that 

/P^(IMNT)  - /tfi^(TNO-)  = tfl^(lMAT)  - ^^(AT)  (B-4) 

I I f I 


and  that 


AH^(2MNT)  - AH^(TNOa)  = /tfI^(2MAT)  - AH^(AT)  (B-5) 

Domalski's  values^  for  ^^(IMAT)  and  AH^(AT)  were  used  in  equation  (4) 
to  give: 

KH  (IMNT)  = ^JI^(IMAT)  - AH^(AT)  + AH^(TNO»)  (B-6) 

f f f 1 " 

= 46.2  - 49.7  + 51  = 48.5 
Rounding  off,  ^SP^CIMNT)  -=  48. 

From  equation  (B-5)  and  Domalski's  values,^ 

^^(2MNT)  = ^^(2MAT)  - ^^(AT)  + AH^(TN02>  (B-7) 

= 50.4  - 49.7  + 51 
= 51.7 

Rounding  off,  AR^(2M)  = 52. 

The  heat  of  formation  of  5-azidotetrazole  (5AZT)  was  obtained  from 
the  assumption  that  AH^(5-azidotetrazole,  solid)  minus  AH^Ctetrazole,  solid) 
is  equal  to  AH^( cyclopentane  azide,  ideal  gas)  minus  (iH^Ccyclopentane, 
ideal  gas) : 

^^(5AZT)  - ^^(CHaN^)  = ^^(CgH^Na)  - (B-8) 

We  inserted  the  values  of  Cox  and  Pilcher*  for  ^^(CgHgNg)  and  /iR^CCgH^p) 

and  Domalski's  value^  for  ^.(CHj,N.)  into  equation  (B-8)  to  obtain 

f 

^^(5AZT)  = AH^(CgHgN3)  - ^^(CgH^g)  + AH^(CHaN4)  (B-9) 

= 52.8  + 18.4  + 56.7  = 127.9 
Rounding  off,  ^^(5AZT)  = 128. 
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Additivity  methods  were  also  used  to  estimate  the  density  of  IMNT 
The  density  is  found  by  dividing  the  molecular  weight  by  the  molar  volume. 
The  molar  volumes  given  in  Table  B-2  show  that  the  difference  between  the 
molar  volumes  of  the  2-isomer  and  the  1-isomer  of  allytetrazole  (in  the 
liquid  state  at  20°C)  is  very  similar  to  the  difference  between  the 
molar  volumes  of  the  2-isomer  and  the  1-isomer  of  vinyl tetrazole  (also 
in  the  liquid  state  at  20°C) . These  values  of  3.4  and  3.5  cra®/mol  are 
within  experimental  error  of  each  other.  Furthermore,  the  slight  change 
in  the  physical  state  for  the  other  pairs  of  Isomers,  namely,  the  liquid 
state  at  25°C  and  the  solid  state,  changes  the  difference  by  less  than 
2%  of  the  average  molar  volumes  to  4.3  and  4,9  cm®/mol.  Therefore  it 
ought  to  be  a very  good  approximation  that  the  difference  between  molar 
volumes  of  the  2-isomer  and  the  1-isomer  for  all  tetrazoles  in  the  solid 
state  is  4,9  cm^/mol. 

From  the  known  molar  volume  of  70.1  cm* /mol  for  2MNT,  the  molar  volume 
of  IMNT  = 70.1  - 4.9  = 65.2  cm*/mol.  The  density  of  IMNT  is  therefore 
115/65.2  =1.76  g/cm*. 
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DISTRIBUTION  LIST 


Quantity  Destlnat Ion 

Sixteen  (16)  Department  of  the  Navy 

Naval  Sea  Systems  Command 
Washington,  D.C.  20362  (N00017) 

Attn;  Code  SEA-0332  (2) 

Code  SEA-0623  (13)  includes  DDC 
Code  SEA-992E  (1) 

One  (1)  Department  of  the  Navy 

Naval  Sea  Systems  Command 
Washington,  D.C.  20362 
Attn:  Code  SEA-0332 

VIA:  Chief,  Defense  Contract  Administration 

Services  Region,  San  Francisco 
866  Malcolm  Road 
Burlingame,  CA  94010  (S0507A) 

One  (1)  Commanding  General 

U.S.  Army  Materiel  Command 
Washington,  D.C.  20315  (DAAGOl) 

Attn:  Library 

Two  (2)  Commanding  Officer 

U.S,  Army  Picatinny  Arsenal 
Dover,  New  Jersey  07801  (DAA21) 

Attn;  SARPA-FR-E-C  (1) 

Technical  Library  (1) 

One  (1)  Commanding  General 

U.S.  Army  Combat  Developments  Command 
Fort  Belvolr,  VA  22060  (DAAK02) 

Three  (3)  Commander,  Naval  Air  Systems  Command 

Washington,  D.C.  20360  (N00019) 

Attn:  Code  AIR-03  (1) 

Code  AIR-350  (1) 

Code  AIR-604  (1) 

Three  (3)  Commander,  Naval  Surface  Weapons  Center/WOL 

Silver  Spring,  MD  20910  (N60921) 

Attn;  Dr.  J.  Enig  (1) 

Dr.  S.  J.  Jacobs  (1) 

Technical  Library  (1) 
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One  (1) 


One  (1) 


Three  (3) 


One  (1) 


Two  (2) 


One  (1) 


commanding  Officer 
Naval  Weapons  Station 
Yorktown,  VA  23491  (N00109) 

Attn:  NEDED 

Commanding  Officer 
Naval  Ordnance  Station 
Indian  Head,  MD  20640  (N00174) 

Attn:  Technical  Library 

Commander 

Naval  Weapons  Center 

China  Lake,  CA  93555  (N60530) 

Attn:  Code  753  (2) 

Code  454  (1) 

Commander 

Naval  Surface  Weapons  Center/OL 
Dahlgren,  VA  22448  (N00178) 

Attn:  Technical  Library 

Air  Force  Systems  Command  (AFSC) 
Andrews  Air  Force  Base 
Washington,  D.C.  20390  (P18600) 

Attn:  SCSM  (1) 

SCT  (1) 

Director 

Advanced  Projects  Research  Agencj' 
Department  of  Defense 
Washington,  D.C.  20301 
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